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not average them when the determined similarity is low, whereby reducing noise without compromise in spatial resolution and in 

temporal resolution. 
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DESCRIPTION 



IMAGE PROCESSING METHOD 

AND IMAGE PROCESSING APPARATUS 

Technical Field 

The present invention relates to an image 
processing method and an image processing apparatus, 
and more particularly to the method and the apparatus 
which can reduce noises existing on an image. 

Background Art 

At present, image processing technology has come 
to be utilized in various fields. 

Image processing is executed in order to cope 
with, e.g., the degradation and amelioration of an 
image acquired by a video tape recorder, a digital 
camera or the like, or it is executed for the purpose 
of, e.g., clearly grasping the pattern or structure 
itself of a structural object in order to inspect if 
the structure is manufactured as designed. 

Also in various medical image diagnostic equipment 
such as an X-ray CT equipment, SPECT equipment and MR I 
equipment, multifarious image processing aspects are 
executed. Regarding the execution of, e.g., the 
depiction of a blood stream or contrast medium stream, 
the sampling of a morbid part, or the contouring of 
internal organs, benefits are extensively accepted. 

The image processing technology consists of 
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various constituent techniques such as a noise 
suppression technique, feature extraction technique and 
pattern recognition technique, and the individual 
techniques are utilized alone or in an appropriate 
5 combination. Incidentally, among such constituent 

techniques, especially a technique which reduces random 
noises contained in an image is an indispensable one 
for more clearly reproducing an object subjected to 
imaging, reconstruction or the like. 
10 Further improvements, however, are required of the 

prior-art image processing technology, especially the 
noise reduction technique. By way of example, so- 
called "smoothing" is widely known as the noise 
reduction technique. The smoothing is such that, when 
15 an input value f (i, j) exists for a certain 

pixel (i, j), an average density in the vicinity of the 
pixel (i, j) is afforded as an output value g(i, j) for 
) the pixel (i, j). Concretely, assuming that n x n 

pixels in the vicinity of the pixel (i, j) are used, the 
20 output value g(i, j) is found as: 

9 (1 ' j) ' J. Jj jb - a ; in5 - c + 1) • f(i + * 3 + 11 

Here, letters a, b, c and d in the above equation (1) 
denote integers. Besides, 1/ (b - a + 1) (d - c + 1) in 
25 the equation (1) is a so-called "weight". 

Incidentally, FIG. 13 shows a case where a, b, c and 
d = -1, 1, -1 and 1. 



Meanwhile, it is generally known that, when the 
average value of n samples independently taken from the 
distribution of a population whose variance is a 2 is 
computed, the variance of the average value becomes 
a 2 /n. According to the equation (1), therefore, the 
"population" and "its variance a 2 " mentioned above 
correspond to a probability distribution whose random 
variable is a component ascribable to noises contained 
in the value of each pixel (i, j), and the variance 
thereof, respectively, so that the component of the 
value f (i, j) of each pixel attributed to the noises 
can be lowered. 

Merely by simply applying the smoothing, however, 
so-called "edge obscurity" appears, the spatial 
resolution of an image is spoilt, and the impression is 
given that the whole image is unsharp. Regarding, for 
example, the medical image referred to above, even in a 
case where a minute blood vessel structure is to be 
depicted with the least possible noise, the noise 
suppression processing based on the equation (1) 
executes the averaging (smoothing) also in pixels which 
do not originally depict the blood vessel structure. 
Therefore, admitting that the noises are suppressed, a 
contrast expressive of the blood vessel structure is 
also lowered by the smoothing, so that the depiction of 
the minute blood vessel structure sometimes becomes 
difficult . 



The present invention has been made in view of the 
above circumstances, and has for its object to provide 
an image processing method and an image processing 
apparatus which can satisfactorily suppress noises 
without giving rise to the unsharpness of an image, and 
which can effectively contribute to other image 
processing techniques, for example, a pattern 
recognition technique . 

Disclosure of Invention 
(Fundamental construction of Present invention) 

The present invention first consists in an image 
processing method characterized in that, in a case 
where a pixel value which is a vector value or a scalar 
value is constructed for a pixel which constitutes a 
certain image, where the goodness of fit between the 
pixel value and another pixel value constructed 
separately from the former pixel value is quantified, 
and where a new pixel value for the pixel is 
constructed by utilizing the other pixel value, the new 
pixel value is constructed by making the contribution 
of the other pixel value large when the goodness of fit 
is high, and by making the contribution of the other 
pixel value small when the goodness of fit is low. 

Besides, in the present invention, the "other 
pixel value constructed separately" may well be 
constructed on the basis of other pixels which 
constitute the image. That is, in this case, the pixel 



values are respectively constructed as "one pixel 
value" and the "other pixel values" each being a vector 
value or a scalar value, for one pixel and the other 
pixels . 

Further, in the present invention, especially in 
case of constructing the new pixel value in such a way 
that a weight function which is a function of the 
goodness of fit is exerted on the goodness-of-f it 
values found for the respective other pixel values, so 
as to determine the weights of the respective other 
pixel values, and that the weighted average of the 
other pixel values is subsequently calculated using the 
weights, the new pixel value is constructed by making 
the weight large when the goodness of fit is high, so 
as to make the contribution of the corresponding other 
pixel value large in the weighted average of the other 
pixel values, and by making the weight small when the 
goodness of fit is low, so as to make the contribution 
of the corresponding other pixel value small in the 
weighted average of the other pixel values. 

According to such processing, the new pixel value 
is constructed with importance attached to the other 
pixel which is decided "similar" to one pixel. 
Therefore, the spatial resolution is not spoilt as is 
in the prior art, or a temporal resolution is not 
spoilt in a case where the image is a dynamic one. 



(Weight function and Weight) 

In the present invention, the weight function 
should preferably be a non-negative monotonously 
increasing function which concerns the goodness of fit. 

Besides, one pixel and the other pixel as 
mentioned above are respectively denoted by x and y, 
and their pixel values by v(x) = (v^ (x) , V2 (x) , 
v K (x)) and v(y) = (y) , v 2 (y) , v K (y)), while the 

new pixel value to be constructed is denoted by v' (x) = 

(v' x (x), v' 2 (x), • . ../ v' K (x)). Further, the goodness 
of fit is denoted by p(x, y) , the weight function is 
denoted by w, and the weight obtained by exerting the 
weight function w on the goodness of fit is denoted by 
w(p(x, y) ) • Then, the processing for taking the 
weighted average should preferably in a form expressed 
as follows: 

]T v k (y) • w (p (x, y) ) 

v k (x) = ^ ' ... (2) 

2^ w (p (x, y) ) 

yeN(x) 

(where N(x) denotes a range in which the other 
pixel is included.) 

(Hereinbelow, the form like the right side of the above 
equation for finding the new pixel value v' (x) shall be 
termed "coherent filter according to the present 
invention" . ) 



(Example of Image processing technology to which 
Present invention is applicable) 

The present invention is capable of various 
applications such as the processing described above, in 
which the new pixel values are constructed over, for 
example, the whole area of the image, thereby to 
reduce the noises of the image or to perform pattern 
recognition in the image. As will be described later, 
the present invention exhibits excellent effects in 
such constituents of the image processing technology. 
(Goodness of fit) 

In the present invention, the goodness of fit may 
well be quantifiable on the basis of a critical 
probability which results from the application of a 
statistical test method to the scalar values respec- 
tively constituting one pixel value and the other pixel 
value stated above. Incidentally, as the "statistical 
testing method" termed here, a "x~square test method", 
for example, can be considered as will be described 
later. Besides, in general, the relationship between 
the critical probability introduced here and the 
goodness of fit shall include especially a case where, 
when one of them increases, the other decreases. That 
is, when the critical probability heightens (the 
goodness of fit lowers) in constructing the new pixel 
value, the weight for the other pixel becomes small 
(the contribution of the other pixel to the new pixel 
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value to be constructed becomes small), and when the 
critical probability lowers (the goodness of fit 
heightens) , the weight for the other pixel becomes 
large (the contribution to the new pixel value to be 
constructed becomes large) . 
(Selection range of "Other pixel") 

Besides, the present invention may well be such 
that the other pixel is selected from within a 
predetermined region around one pixel mentioned above, 
and especially that the predetermined region is the 
vicinity centering round one pixel mentioned above. 
With such a construction, the other pixel decided 
similar to one pixel is usually much likely to exist 
around one pixel mentioned above, so that, when the 
region as stated above is selected, wasteful 
calculations can be omitted in effectively finding the 
new pixel value. 

(Method of Constructing Pixel values) 

In the present invention, various methods of 
constructing the pixel values (both one pixel value and 
the other pixel value) can be selected. More 
concretely, in the first place, in a case where the 
image consists of a plurality of images, the pixel 
values can be set as vector values in each of which the 
scalar values of pixels at the same points in all of 
the plurality of images are arrayed. 

On this occasion, the plurality of images may well 



be a plurality of static images which constitute a 
dynamic image. In such a case, the present invention 
does not spoil the spatial resolution and does not 
spoil the temporal resolution, either, in finding the 
new pixel value. Incidentally, a favorable application 
in such a case is, for example, a case where the 
dynamic image is a dynamic CT image acquired by a 
medical image diagnostic equipment. 

Besides, the second method of constructing the 
pixel values corresponds to a case where the image 
consists of a plurality of sorts of images which 
concern an identical subject. Herein, as in the 
foregoing, the pixel values can be set as vector values 
in each of which the scalar values of pixels at the 
same points in all the plurality of sorts of images are 
arrayed. Incidentally, a favorable application in such 
a case is, for example, a case where the image is an 
image acquired by multiwindow radiography in a SPECT 
equipment or PET equipment, and where the plurality 
of sorts of images are images different from one 
another, based on gamma rays emitted from different 
radioisotopes. Another favorite application is, for 
example, a case where the image is a color image, and 
where the plurality of sorts of images are images 
different from one another as acquired by decomposing 
light into three primary colors. 

Further, as the third method of constructing the 



pixel values, the image is a single image, and the 
pixel values can be constructed on the basis of the 
single image. More concretely, the pixel values can be 
constructed as, for example, vector values in which the 
scalar values of a certain pixel on the single image 
and the scalar values of a pixel existing in the 
vicinity of the certain pixel are respectively arrayed. 
Incidentally, according to such a method of 
constructing the pixel values, the present invention is 
applicable to any digital image. 

(Preferred apparatus construction for Realizing Image 
processing of Present invention) 

Lastly, an apparatus construction for realizing 
the image processing method described above should 
preferably comprise goodness of fit quantification 
means for deciding the goodness of fit between a pixel 
value constructed as a vector value or a scalar value 
and another pixel value separately constructed, as to a 
pixel which constitutes a certain image, and pixel 
value calculation means operating in constructing a new 
pixel value concerning the pixel by utilizing the other 
pixel value, for constructing the new pixel value by 
making the contribution of the other pixel value large 
when the goodness of fit is high and by making the 
contribution of the other pixel value small when the 
goodness of fit is low. 
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Brief Description of Drawings 
FIG. 1 is a table showing the consequent 
relationships of a critical probability p(x, y) , a 
weight w(p(x, y) ) and a null hypothesis. 
5 FIG. 2 is a table showing various examples of 

applications in embodiments of a coherent filter 
according to the present invention, by classifying them 
in accordance with differences in a method of 
constructing a pixel value v(x). 
10 FIG. 3 is a diagram showing an example of the 

schematic arrangement of an X-ray CT equipment in which 
the coherent filter is incorporated. 

FIG. 4 is an explanatory diagram for conceptually 
explaining image processing based on the coherent 
15 filter according to the present invention. 

FIG. 5 shows an example in which a plurality of 
images constituting a dynamic CT image are subjected to 
} noise suppression processing based on the coherent 

filter according to the present invention; wherein (a) 
20 shows an original image not processed yet, and (b) and 

(c) show processed images, respectively. 

FIG. 6 is a flow chart showing the flow of the 
noise suppression processing based on the coherent 
filter in an embodiment. 
25 FIG. 7 shows an example in which a single CT image 

is subjected to noise suppression processing based on 
the coherent filter according to the present invention; 



wherein (a) shows an original image not processed yet, 
(b) shows a processed image, and (c) shows an image 
obtained by averaging the image shown in (a) and the 
image shown in (b) . 

FIGS. 8(a), (b) , (c) and (d) are enhanced images 
each showing an example of 3D sampling of blood 
vessels, and FIG. 8(e) is an image obtained by sampling 
the gray matter of a brain as another sampling example. 

FIG. 9 is an explanatory diagram showing one 
example of a method of constructing a pixel value v(x) 
from a single image. 

FIG. 10 is an explanatory diagram showing another 
example of the method of constructing a pixel value 
v(x) from a single image. 

FIG. 11 is a graph showing a time-density curve in 
the case of combining a coherent regression method and 
an AUC method, and the curve in the case of combining a 
prior-art method and the AUC method. 

FIG. 12 relates to another estimation method for a 
variance a 2 in Equation (4), etc. in the body of the 
specification, and is an explanatory diagram for 
explaining a calculation method for a variance a 2 + + 
based on the other estimation method. 

FIG. 13 is an explanatory diagram for explaining 
"'smoothing processing" which is a noise suppression 
method in the prior art. 



Best Mode for Carrying Out the Invention 
Now, embodiments of the present invention will be 
described with reference to the drawings. 
Incidentally, hereinafter, more general description 
(Item No. "I" below) on the outline of a "coherent 
filter" according to the present invention shall be 
first made, and examples (Item Nos. "II" and "III" - 
"VIII") in which the coherent filter is applied to 
various image processing aspects shall be thereafter 
described in succession. 

The description of such various examples of 
applications shall be further followed by description 
on the most general form of the coherent filter in 
the embodiment and an example of application (Item 
No. "IX") thereof, and the supplementary items (Item 
No. "X") of the embodiments shall be lastly summed up. 
(I General description of Coherent filter according to 
Present invention) 

First, there will be made preparatory description 
for describing the "coherent filter" according to the 
present invention. 
(1-1 Pixel value v(x)) 

In general, a digital image acquired through 
imaging means such as a camera is constituted by a 
plurality of pixels (or, the image can be considered a 
set of such pixels) . In the ensuing description, the 
position of each of the pixels will be expressed as a 



vector x (that is, a vector of coordinate values), a 
value (for example, a numerical value expressive of a 
density) which the pixel x has will be expressed as a 
K-dimensional vector. In case of a two-dimensional 
image, the pixel x is a two-dimensional vector which 
indicates the coordinate values (x, y) expressive of 
the position on the image. A "pixel value v(x)" 
defined for a certain pixel x is denoted as: 

v(x) = (v 1 (x) / v 2 (x), v K (x)) ... (3) 

Incidentally, v 1 (x) , v 2 (x) , ••• and v K (x) on the right 
side of this equation (3) shall be respectively called 
"scalar values" concerning the pixel x stated below. 

By way of example, when the image is a "color 
image' 7 , each pixel has brightnesses (scalar values) in 
three primary colors (red, green and blue), so that the 
pixel value v(x) of each pixel can be considered as a 
vector whose number of dimensions is K = 3 (suppose a 
case where, in the right-side terms of the above 
equation (3), suffixes are, for example, "red' 7 , "green" 
and "blue", and refer to Equation (3. 1) to be 
mentioned later) . Besides, in a case, for example, 
where the image is a dynamic image consisting of K 
static images, and where each pixel of the nth image 
has a scalar value v n (x), a K-dimensional vector 
value v n (x) = (v x (x) , v 2 (x) , • --, v K (x)) which is 
constructed by arraying the pixel values (scalar 
values) of pixels x at the same common points (same 
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coordinates) on the K static images is the pixel value 
being the vector value to be described later. 
(1-2 Goodness of fit or Critical probability p(x, y) 
and Weight w(p(x, y) ) ) 
5 For the pixel x, a suitable set N(x) of pixels is 

considered (the set N(x) may well include the pixel x) . 
Subsequently, a weight w(p(x, y) ) is considered between 
the pixel x and each pixel y which is the element of 
the set N(x) . The weight w(p(x, y) ) has properties as 
10 stated below. 

(1-2-1 Goodness of fit or Critical probability 
P(x, y) ) 

First, there will be explained the significance of 
the function p (x, y) which governs the value of the 

15 weight w(p(x, y) ) . This function p(x, y) is means for 

quantifying the "goodness of fit" as termed in the 
present invention. Generally speaking, it gives a 
} concrete numerical value indicating the degree to which 

the pixel x and a pixel yeN(x) are similar in any sense 

20 (for example, the degree of a statistical difference 

which is noticed between the pixel values v(x) and v(y) 
of both the pixels x and y) . 

More concretely, by way of example, when the 
function p (x, y) gives a small value, it is judged that 

25 a "statistically significant difference does not exist 

(= the goodness of fit is high)" between the pixel 
values v(x) and v(y), and that the pixels x and y are 



much likely to be similar, and when the function 
p(x, y) gives a large value, it is judged that a 
"statistically significant difference exists = the 
goodness of fit is low)". 

Meanwhile, it must be considered that noises are 
inevitably contained in the pixel values v(x) and v(y) 
(or the scalar values vi (x) , v K (x) and vi (y) , 

v K (y)). By way of example, when a case is considered 
where the image has been acquired by a CCD imaging 
device, noises which are ascribable to dark currents 
within the device and the irregular fluctuations of the 
quantity of light entering from the outside exist as to 
the individual pixels constituting the image. 

In general, such noises take diverse values as to 
all the pixels. Therefore, even if the pixels x and y 
reflect an identical object (in the outside), they 
sometimes fail to have the same values on the image 
which is actually observed. Conversely speaking, 
supposing a situation where the noises have been 
removed from the pixels x and y both of which reflect 
the identical object, they are displayed (= recognized) 
on the image as representing the identical object, and 
both the pixels ought to have the same (or very close) 
pixel values. 

Therefore, when the concept of a "null hypothesis" 
well known in a statistical testing method is employed 
as to the function p(x, y) in consideration of the 



properties of the noises stated above, this function 
p(x, y) can be concretely elucidated as below. In case 
of building up the null hypothesis H that "the pixels x 
and y have the same pixel values when their noises have 
been removed", in other words, that "v(x) = v(y) except 
differences ascribable to the noises of both the 
pixels'' (that is, when such a proposition is concluded, 
it is considered that "both the pixels x and y are 
similar (= the goodness of fit is high) ", the function 
p(x, y) can be said a critical probability (or a 
significant level) in the case where the hypothesis H 
is abandoned (in this case, p(x, y) is defined as the 
function whose range is [0, 1] (p(x, y)e[0, 1])). 

Accordingly, in a case where the critical 
probability p(x, y) is high, that is, where the 
abandonment is much likely to be erroneous, the 
possibility of fulfilling the hypothesis H can be said 
high, and conversely, in a case where the critical 
probability p(x, y) is low, that is, where the 
abandonment is little likely to be erroneous, the 
possibility of failing to fulfill the hypothesis H can 
be said high. (Incidentally, as is well known in the 
statistical testing, even when the hypothesis H is not 
"abandoned", it is not meant that this hypothesis is 
"true". In this case, it is merely meant that the 
proposition indicated by the hypothesis H cannot be 
negated. ) 



18 

(1-2-2 Weight w(p(x, y) ) ) 

Meanwhile, as apparent from its expression, the 
weight w(p(x, y) ) is the function of the critical 
probability p(x, y) as explained above (more generally, 
it can be constructed as w(p(x, y) ) (where p(x, y) 
denotes the goodness of fit) as the function of the 
goodness of fit) . Besides, the weight function w which 
is exerted on the critical probability p (x, y) obtained 
for each combination of the pixels x and y, in order to 
find the weight w(p(x, y) ) has, in general, the action 
of incarnating the above ^abandonment". Concretely, 
the weight function w is adjusted in such a manner 
that, in a case where the critical probability p (x, y) 
is high, the value of the weight function w, namely, 
the weight w(p(x, y) ) takes a large positive value, and 
that, in the converse case, it takes a small positive 
value (or "0") (the concrete form of the weight 
function w will be explained later) . That is, the 
weight w(p(x, y) ) takes the large value in a case 
where the pixels x and y are likely to fulfill the 
proposition indicated in the hypothesis H, and it takes 
the small value in the converse case. By way of 
example, the function w may well be especially 
constructed so as to be capable of taking only the two 
values of "0" and a constant value which is not "0". 

Incidentally, the relationships among the 
hypothesis H, critical probability p(x, y) and weight 
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w(p(x, y) ) thus far explained are summed up in FIG. 1. 
Besides, the weight function w(t) can be more generally 
said a "non-negative monotonously increasing function 
defined for te[0, 1]", and a property which this 
5 function w(t) is to satisfy may be, at least, such. 

(1-3 Coherent filter) 

Owing to the preparatory description thus far 
made, the "coherent filter" according to the present 
invention is derived as stated below. First, as to a 

10 certain pixel x constituting an image, the weights 

w(p(x, y) ) stated above are computed for all the pixels 
y being the elements of a set N(x) . Subsequently, 
using the plurality of weights w(p(x, y) ) , new scalar 
values v' k (x) to constitute the pixel x are computed by 

15 the following equation: 

2 v k(Y)w(p(x, y) ) 
, , . yeN(x) 

v k (x) = ^ v ' ■ - - (4) 

2^ w(p(x, y) ) 

yeN(x) 

Here, k = 1, 2, K. Besides, using v'^fx) 

obtained with this equation, a pixel value (new pixel 
20 value) v' (x) after the transformation of the pixel x is 

constructed as: 

v' (x) = (v'!(x), v' 2 (x), v' K (x)) (5) 

Herein, a filter which transforms the pixel value 
v(y) = (v 1 (y) / v 2 (y), v K (y)) (including a case 

25 where y = x holds) as expressed by the above equation 

(4), into V (x) = (v'ifx), v' 2 (x), v' K (x)), is the 



form of the "coherent filter" according to the present 
invention, in this embodiment. As apparent from the 
expression, the new pixel value represents the weighted 
average values of the scalar values v^ty) constituting 
the pixel value. 

Such processing brings forth results as stated 
below. The pixel value v' (x) comes to express a vector 
which is constituted by the weighted average values 
v') c (x) with importance attached to the pixel y that may 
probably take the same pixel value as that of the pixel 
x except noises (= that is much likely to fulfill the 
proposition of the hypothesis H) . Besides, if such 
pixels y exist in a sufficient number, the pixel value 
v' (x) comes to have a value with only the noises 
suppressed by the above averaging action, without 
deviating from the true value of the pixel x as ought 
to be had. 

By the way, even in a case where the critical 
probability p(x, y) is low, and the null hypothesis H 
is accordingly "abandoned", so that the weight 
w(p(x, y) ) becomes small, the hypothesis is not always 
completely "abandoned" as understood from the above 
description or the expression of Equation (3) . Even in 
a case where the critical probability p(x, y) is close 
to "0" ( = 0%), the weight w(p(x, y) ) * 0 may well be 
set (at a positive value which is smaller than in a 
case where the probability p(x, y) is close to "1") , 



though such depends upon the concrete form of the 
weight function w to be described later. 

(Incidentally, p(x, y) = 1 corresponds to a case where 
v(x). = v(y) holds as will be stated later.) That is, 
unlike the complete abandonment, a small contribution 
may well be admitted. (By the way, if the weight 
w(p(x, y) ) = 0 is set in such a case, this has the same 
meaning as the complete abandonment. Refer to 
Equation ( 14 ) later . ) 

Such processing can be generally elucidated as 
below. It can be said an image processing method in 
which, when (a plurality of) pixels x constituting a 
certain image exist, the goodness of fit between each 
pixel x and a certain desired pixel y (set as yeN(x) in 
the foregoing) is quantified (on the basis of p(x, y) 
in the foregoing) , and in which, when the goodness of 
fit is high, a large contribution is admitted for the 
pixel y in the weighted averaging process utilizing the 
pixel value v(y), and when the goodness of fit is low, 
only a small contribution is admitted, whereby the 
noises of the pixel x are effectively suppressed. In 
other words, when the pixels x and y are "similar to 
each other", the pixel y is caused to contribute to the 
averaging process more, and when they are "dissimilar 
to each other", the pixel y is almost or utterly 
ignored. 

By executing such processing over the whole image, 
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a very high effect of suppressing noises can be 
exhibited without almost no unsharpness of the image. 
Moreover, not only in the application to the noise 
suppression, but also in the field of, for example, 
5 pattern recognition, excellent effects can be exhibited 

by bringing the weight function or the coherent filter 
into a preferred concrete form. Incidentally, how 
these effects are in the concrete will be stated in the 
ensuing description on various image processing 
10 aspects. 

(II Various examples of Applications; Classification 
based on Differences in Method of Constructing Pixel 
value v (x) ) 

Meanwhile, the present invention is applicable to 
15 multifarious fields on the basis of the general form of 

the coherent filter described above. The applications 
can be summed up as listed in FIG. 2 by way of example. 
) By the way, in the list of FIG. 2, the various examples 

of applications in multifarious fields are classified 
20 in accordance with differences in the method of 

constructing the pixel value v(x). According to 
FIG. 2, it is understood that the pixel value v(x) 
being the vector value stated above can be constructed 
by individual schemes "from a plurality of images" 
25 (obtained by, for example, one time of radiographing) , 

"from a plurality of sorts of images' 7 or "from a single 
image", and that the various examples of applications 
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concerning this embodiment are classified into the 
individual schemes . 

Now, there will be successively described the 
various examples of applications in which the coherent 
5 filter having the general form as stated above is 

applied in more practicable scenes, along the 
classification table shown in FIG. 2, and how three 
sorts of methods of constructing the pixel value v(x) 
shown in FIG. 2 are realized in the concrete (in an 

10 different from that shown in FIG. 2) . 

(Ill Case 1 of Applying Present invention to X-ray CT 
radiography; Dynamic CT radiography (To eliminate 
Noises from Plurality of images) ) 

First, there will be described a case where the 

15 coherent filter is applied to so-called "dynamic CT" 

radiography in an X-ray CT equipment. Here, . as shown 
in FIG. 3, the X-ray CT equipment 100 includes an X-ray 
) tube 101 as well as an X-ray detector 102, a data 

acquisition system 103 (DAS) , a preprocessing unit 104, 

20 a memory unit 105, a reconstruction unit 106 and an 

image display unit 107. Also included are a control 
unit 108 which controls these constituents, and an 
input unit 109 which is capable of various settings and 
inputs of an X-irradiation condition, a radiographing 

25 mode, etc. 

According to the X-ray CT equipment 100, a subject 
or patient P is exposed to X-rays while the X-ray tube 



24 

101 and the X-ray detector 102 are being rotated around 
the patient P, whereby the internal state of the 
patient P can be observed as a CT image, for example, 
tomogram. 

5 On this occasion, the constituents mentioned above 

operate as outlined below. X-rays emitted from the 
X-ray tube 101 by the voltage application of a high 
voltage generator 101a and transmitted through the 
patient P, are converted into an analog electric signal 

10 by the X-ray detector 102. Thereafter, the analog 

electric signal is subjected to a digital conversion 
process by the data acquisition system 103, and various 
correction processes by the preprocessing unit 104, and 
the resulting digital signal is stored in the memory 

15 unit 105 as projection data. Besides, the tomogram or 

any other CT image is reconstructed by the reconstruc- 
tion unit 106 on the basis of the projection data 
) obtained in such a way that the X-ray tube 101 and 

the X-ray detector 102 has made, for example, one 

2 0 revolution (= 360°) round the patient P, and the image 

display unit 107 displays the CT image. Incidentally, 
it is also possible to store the reconstructed CT image 
in a storage device 10M. 

Here, the "dynamic CT" radiography mentioned above 

25 signifies a radiographing scheme wherein the X-ray tube 

101 and the X-ray detector 102 iteratively radiograph 
the identical part of the patient P (by iterative 
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scans, and iterative radiographing based on continuous 
rotation is often performed in a continuous rotating CT 
equipment) so as to successively acquire projection 
data, and reconstruction processes are successively 
performed on the basis of the projection data so as to 
obtain a series of images in time series. (In this 
case, an image display in the image display unit 107 is 
controlled by, for example, an unshown counter so as to 
be presented a predetermined time period after a scan 
start point or end point concerning the acquisition of 
the source projection data of the pertinent image.) 

Accordingly, the image thus acquired and displayed 
becomes a so-called "dynamic image" which consists of 
a plurality of time-series static images as in a 
motion picture or the like. Incidentally, such a 
radiographing scheme is typically used in order to 
analyze the morbidity of a morbid part, for example, 
constriction or obstruction in a blood vessel by 
injecting a contrast medium into the patient P and 
observing and investigating the change of the image 
with the passage of time. Besides, a scheme wherein 
the CT radiographing of an identical part is performed 
only twice before and after the contrast medium is 
given can be considered as dynamic CT radiography in a 
broad sense. 

Meanwhile, in the prior art, in a case where the 
patient P has undergone any change (for example, the 
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density change of the contrast medium, or a respiratory 
motion is generally considered) in the "dynamic CT" 
radiographing mode as stated above, for example, while 
K times of radiographing are performed, smoothing in 
5 the direction of time has been inevitably executed in 

order to suppress image noises without spoiling a 
spatial resolution- Consequently, it has been 
impossible to avoid the drawback that a temporal 
resolution is spoilt. 

10 In this regard, the image acquired by the dynamic 

CT radiography is the dynamic image and is intended to 
minutely observe temporal changes, as stated above, so 
that the spoilage of the temporal resolution cannot be 
originally said a favorable situation. 

15 When the coherent filter according to the present 

invention is utilized, it is possible to execute 
processing to be stated below (hereinafter, called 
) "dynamic coherent filtering'') which can suppress noises 

for all of K static images (a plurality of images) 

20 without spoiling a temporal resolution. 

First, regarding a pixel x which is defined as to 
the K static images constituting the dynamic image 
obtained in the above way, the following can be 
constructed as a pixel value v(x) as already stated: 

25 v(x) = (vx(x), v 2 (x), v K (x)) 

• • • (3 mentioned before) 
Here, suffixes 1, 2, •*•, and K at the respective terms 
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on the right side are Nos . allocated through the K 
static images (refer to the foregoing description on 
Equation (3) ) . 

Subsequently, the concrete form of a weight 
function wl in this case is given by, for example, the 
following equation : 



wl(p(x, y) ) = exp 



y (vk(x) ~ Vk(Y) y 
k ~ K(2a k )2 



(6) 



Here, yeN(x) holds, and the set N(x) may be set at 

10 will as to the pixel x (= it may be set in accordance 

with any criterion) . In actuality, however, it can be 
said that the pixel x and the pixel y lying at a 
position far distant from the pixel x fulfill the 
hypothesis "v(x) = v(y) except differences ascribable 

15 to the noises of both the pixels" generally at a 

low possibility, and hence, it has a practical 
significance, such as enhancement in a computing speed, 
that the set N(x) is limited by the criterion of being 
a set of pixels close to the pixel x. 

20 Accordingly, as an example here, the set N(x) 

is assumed to be a set of pixels which are included 
in a rectangular area centering around the pixel x. 
More concretely, in a case, for example, where all 
pixels constituting one static image now noticed are 

25 128 x 128 pixels, the set N(x) may be selected as the 

area of 3 x 3 pixels around the pixel x. Also, in a 
case where all pixels are 512 x 512 pixels, the set 
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N(x) may be selected as the area of 13 x 13 pixels 
around the pixel x. 

Besides, in the above equation (6) denotes the 
standard deviation of noises estimated assuming that 
the individual pixels of the kth static image have the 
noises at a certain common degree, while C denotes a 
parameter for permitting to determine and adjust the 
degree of an action in the case where the weight 
wl(p(x, y) ) is substituted into Equation (4) mentioned 
before . 

Now, these quantities and C will be explained 
in succession. 

First, the quantity in Equation (6) will be 
explained (it will be explained as a variance a^ 2 
below) . As stated above, the variance a^ 2 is that of 
the noise components which the scalar values of the 
individual pixels on the kth static image have. 
Besides, the variance a^ 2 in the above equation (6) has 
been estimated assuming that the scalar values of the 
individual pixels of the kth image contain the noises 
of the variance a^ 2 being a constant value. In 
general, such an assumption is satisfactorily 
reasonable with a background as stated below. 

First, in a state where the size of the patient P, 
and the structures of the X-ray tube 101 as well as the 
X-ray detector 102, the reconstruction unit 106, etc. 
are constant, and where the energy of irradiating 
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X-rays is held constant, the noises of a CT image are 
determined by the quantity of X-irradiation, that is, 
the product between a tube current in the X-ray tube 
101 and an irradiating time period (what is called 
5 "current time product'' (mA-s) ) as is in a proportional 

relation with the quantity of X-irradiat ion . On the 
other hand, it is also known that the noises of the CT 
image are additive and substantially conform to a 
Gaussian distribution. More specifically, when v n °(x) 

10 is let denote a true value (a value obtained by 

eliminating the contribution component of noises) as to 
any scalar value v n (x) (n = 1, 2, K) which 

constitutes the pixel value v(x) of a certain pixel x, 
the value of the difference between them; v n (x) - 

15 v n ° (x) substantially conforms to a Gaussian 

distribution of average 0 and variance a^ 2 
(incidentally, the quantity of X-irradiation or the 
) current time product mA-s and the variance a^ 2 of the 

noises are substantially in an inversely proportional 

20 relation) . 

Besides, although the variance a^ 2 depends also 
upon the position itself of the pixel x (for example, 
coordinate values x = (x, y) as stated above) , it is 
negligible in the ordinary X-ray CT equipment 100 

25 because a physical X-ray filter for adjusting the 

quantity of X-irradiation (what is called "wedge" or 
"X-ray filter" constructed of, for example, copper foil 
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or metal lump) (not shown) is interposed between the 
X-ray tube 101 and the X-ray detector 102. The reason 
therefor is as stated below. The wedge has the 
function of adjusting (absorbing or shielding) part or 
all of the irradiating X-rays so that approximately the 
same quantity of X-irradiation may be detected in any 
X-ray detector 102, by utilizing the fact that the 
patient P is constructed of substances being substan- 
tially equal in density to water. Accordingly, such a 
wedge results in producing the effect of making the 
variance a^ 2 of the noises a constant value which 
hardly depends upon the position of the pixel x. 
(Incidentally, the wedge is disposed usually with the 
original object of effectively utilizing the dynamic 
range of the X-ray detector 102.) 

In the light of the above facts, it is reasonable 
that, in the K static images acquired by the dynamic CT 
radiography, the variance a^ 2 being the constant value 
is estimated for all the pixels on the kth static 
image. Of course, the expansion of this embodiment 
will be easily suggestible as to a case where the 
variance differs every pixel (as will be described in 
detail in Section X-3 later) . 

Subsequently, it becomes a problem what numerical 
value is to be used as the variance o^ 2 in concretely 
calculating the above equation (5) . Such a problem is 
ascribable to the fact that, even if the shape of the 
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distribution of noises (in the above, the Gaussian 
distribution) can be ordinarily supposed, the concrete 
value of the variance a^ 2 is often unclear. 

Further, in general, the radiography may well be 
5 performed by altering the quantity of irradiation 

(X-ray tube current x irradiating time period (mA-s)) 
every time of radiographing. 

Now, letting a^ 2 denote the variance of noises 
which the scalar values of individual pixels have in 

10 the kth image (k = 1, 2, K) , and denote the 

quantity of irradiation used in the radiographing of 
the kth image, a^ 2 is proportional to R^. Accordingly, 
when akg 2 can be designated for at least one image k = 
kg, the variance a^ 2 can be precisely estimated 

15 concerning also the other images k's, in accordance 

with: 

R k0 

In this embodiment (to which such circumstances 
20 apply) , the concrete numerical value of the variance 

a^ 2 can be estimated by a method stated below, as to at 
least one image k. 

Effective is a method wherein, using N (1 < N < K) 
images as to which it can be assumed that the patient P 
25 hardly underwent changes in K times of radiographing, 

an expected value Eta^ 2 ] for the variance o^ 2 is found 
by an actual measurement. For the brevity of the 
ensuing description, it is assumed that the N images 
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are based on the same quantities of irradiation, so the 
variance a^ 2 is constant (written as a 2 ) for k = 1, 
2, •••/ N. It is conjectured that noises which are 
contained in individual scalar values v^ (Xf ) , 
V2 (xf ) , and v K (xf) constituting the pixel value 

v(xf) of a certain pixel Xf in the N images conform to 
the Gaussian distribution of average 0 and variance a 2 
as stated before. Therefore, using the average value 
of the scalar values as follows: 

1 N 

v *< x f> = - 5> k (x f ) ■•• (7) 

k = l 

the expected value E[a 2 ] for the true variance a 2 can 
be calculated as: 

N 2 
15 E[a 2 ] = — ±— J) {v k (x f ) - v*(x f )} • ■ • (8) 

N " 1 k = l 

Herein, this expected value E[a 2 ] of the variance can 
J be considered reasonable for all the pixels x on all 

the K static images as stated before, and it is a value 
as to which a certainty is guaranteed at or above a 
20 predetermined degree in the use thereof as a substitute 

for the true variance a 2 . In the actual calculation 
of the above equation (6), accordingly, the expected 
value E[a 2 ] may be substituted into the term a 2 of 
Equation { 6) . 

25 More concretely, such an expected value E[a 2 ] may 

well be calculated from actual measurement values based 
on, for example, the first and second static images 
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among the K static images (in the above equations (7) 
and (8"), this corresponds to setting N = 2). Besides, 
regarding the pixel Xf for use in the actual 
calculations of Equations (7) and (8), such a 
5 contrivance may well be made that only the suitable 

pixel Xf except parts where the air or bones, for 
example, are imaged is selected (in a case where a 
plurality of pixels have been selected, all the 
expected values E[a 2 ] obtained are averaged). Further, 

10 in general, it is still better to make such a 

contrivance as suppressing influence ascribable to the 
movement of the patient P. (Incidentally, the 
evaluation of the variance a 2 of noises in Equation (6) 
will be referred to again in "Supplementary items of 

15 Embodiments' 7 to be stated last. 

It will be easily suggestible that, even in a case 
where the quantities of irradiation are not constant 
} in the radiographing steps of the N images, the 

variance Op 2 is precisely estimated by utilizing 

20 the proportionality thereof to the quantity of 

irradiation R^. 

Next, the parameter C in the above equation (6) 
will be explained. First, in Equation (6), the idea of 
the critical probability p(x, y) stated concerning the 

25 general form is included as below. An expression 

within a root in the numerator of the right side of 
Equation (6) agrees with a x 2 value which conforms to a 
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so-called "x~ sc 3 u are distribution" . A value which is 
obtained in such a way that the %^ value is divided by 
(2a) 2 and that the entirety in braces is put at the 
shoulder of e, is the critical probability pi (x, y) 
itself. That is: 



pl(x, y) = A exp 



Besides, the above equation (6) is nothing but the 
following formed as to pi (x, y) expressed by this 
10 equation (9) : 

_C 

wl(p(x, y) ) = {pl(x, y)} 2 • • - (10) 

Incidentally, letter A denotes a constant which 
normalizes the critical probability pi so as to become 
a value of (0 through 1) . 

15 After all, the critical probability p(x, y) stated 

above as to the general form is not positively 
indicated in Equation (6), but the actual state of 
the weight wl(p(x, y) ) can be deemed the very 
function (Equation (10) ) of the critical probability 

20 (= pi (x, y) ) as described above, that is, it is the 

''function of the goodness of fit". (As stated before, 
the critical probability and the goodness of fit are in 
the relationship that, when one increases, the other 
increases, too.) 

25 Besides, as seen from the above equation (10) , 

the parameter C has the effect of determining how 
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sensitively the weight wl(p(x, y) ) reacts with the 
critical probability pi (x, y) . That is, when the 
parameter C is enlarged, the weight wl(p(x, y) ) 
approaches zero merely by the slight decrease of the 
5 critical probability pi (x, y) . In contrast, when the 

parameter C is made small, such an oversensitive 
reaction can be suppressed. Incidentally, the 
parameter C may be concretely set at about 1 through 
10, and it should preferably be set at C = 3 . 

10 In this embodiment, the similarity decision 

concerning both the pixels x and y, in other words, the 
decision of the abandonment of the foregoing null 
hypothesis H concerning both the pixels x and y, is 
made by a so-called "x~ sc i uare test method" (statistical 

15 testing method) on the basis of the critical 

probability pi (x, y) as understood from the above 
description . 

Besides, as seen from the expression of the above 
equation (6), in the present invention, the steps of 

20 computing the critical probabilities p(x, y) for the 

respective combinations of the pixels x, y and 
thereafter finding the weights w(p(x, y) ) need not 
always be traced, but a composite function (wOp) may 
well be directly computed without concretely 

25 calculating the critical probabilities p(x, y) . 

As thus far described, the variance a 2 is 
estimated (for example, E[o 2 ] in Equation (8)), and the 



parameter C is appropriately determined (for example, 
C = 3), whereby the concrete weights wl(p(x, y) ) can be 
found for all the pixels y which are included in the 
set N(x) defined as to the certain pixel x (for 
example, the area corresponding to 3 x 3 pixels around 
the pixel x as stated above), by using Equation (6). 
Thereafter, the weights wl(p(x, y) ) are used instead of 
w(p(x, y) ) in the above equation (4), whereby the 
concrete numerical calculations of the coherent filter 
can be executed. As a result, the pixel value v' (x) = 
(v' i (x) , v'2< x )/ •"*/ v' K (x)) (= Equation (5)) with the 
noises intensely suppressed, in other words, such K 
static images or such a dynamic image can be obtained 
without spoiling, not only the temporal resolution, but 
also the spatial resolution. 

It is FIG. 4 that shows such image processing so 
as to be easy of being conceptually grasped. First, 
referring to FIG. 4 (a) , in first, second, ... and 
Kth static images, a rectangular area N 3x3 (x) 
centering around a certain pixel x and corresponding to 
3x3 pixels is supposed for the pixel x. Letting 
denote the pixel at the left upper corner of the 
rectangular area N 3x3 (x), this pixel has a pixel 
value v(y^) as also indicated in FIG. 4. 

Herein, a weight wl (p (x, y^) ) is computed by the 
above equation (6) on the basis of scalar values 
v l(Yl)/ v 2(¥l)' '•"/ and v K(yi) constituting the pixel 
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value v(yi), and scalar values (x) , V2 (x) , and 
v K (x) in a pixel value v(x) { FIG . 4(b)), Besides, the 
remaining pixels Y2f *' # / and Y8 of the rectangular 
area N 3x3 (x) are similarly processed, so that weights 
5 wl(p(x, yi)), •••/ wl(p(x, y 8 ) ) andwl(p(x, x) ) are 

finally obtained as shown in FIG. 4(b). (In this case, 
a critical probability p(x, x) is xx l" in accordance 
with Equation (9), so that the weight wl (p (x, x) ) is 
also "1" in accordance with Equation (10) (= the 
10 maximum weighting is done) .) 

Subsequently, the weights wl(p(x, yi)), 
wl(p(x, ys)) and wl(p(x, x) ) thus obtained are 
respectively multiplied by the scalar values V] c (y^), 

v k(Y2)' '""/ v k(Y8)' and v k( x ) of the corresponding 
15 pixels in the kth image, the resulting products have 

their summation taken (corresponding to a numerator in 
the above equation (4)), and the resulting summation is 
/ divided by the summation of the weights wl concerning 

the rectangular area N 3x3 (x) (corresponding to a 
20 denominator in the same equation (4)). Then, a scalar 

value v'j c (x) with noises suppressed, concerning the 
pixel x in the kth image can be calculated (FIG. 4(c)). 
Besides, scalar values v'^fx) with noises suppressed 
are calculated for all the images of k = 1, 2, 
25 K by using the same weights wl(p(x, yi)), 

wl(p(x, ys)) and wl(p(x, x) ) , whereby pixel values 
v'fc(x) = (v'i(x), v f 2( x )/ • • • , v'k(x) ) with noises in. 
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the pixel x suppressed are obtained. When the above 
calculations are iterated for all the pixels x, K 
images with noises suppressed are obtained. 

Images which are constituted by the pixel values 
5 v' (x) thus calculated become, for example, as shown in 

FIG. 5. FIG. 5(a) shows an original image (k = 11) 
before the application of the coherent filter, among 
K = 23 static images in all, and (b) and (c) show 
images after the application of the coherent filter, 

10 for k = 11 and 23, respectively. Random noises seen in 

the original image of FIG. 5(a) are sufficiently 
suppressed in the images of FIGS. 5(b) and (c) . 

Incidentally, the processing steps thus far 
explained may be executed in accordance with, for 

15 example, flow charts as shown in FIG. 6. Besides, in 

order to realize the calculations, image display, etc. 
according to the processing steps, on the actual X-ray 
) CT equipment 100, they may be executed by disposing, 

for example, an image processing unit 110 which is 

20 configured of a variance estimation portion 111, a 

weight calculation portion 112 and a pixel value 
calculation portion 113 as shown in FIG. 3. 

Among them, the weight calculation portion 112 is 
so constructed that the weight wl(p(x, y) ) is directly 

25 calculated from the pixel values v(x) and v(y) as in 

the foregoing procedure. Accordingly, the calculation 
portion 112 is a device which calculates the weight 
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directly without concretely calculating the value of 
the critical probability pi (x, y) (that is, by 
including a critical probability calculation portion 
(a "goodness of fit quantification portion" termed in 
5 the present invention) ) . Incidentally, unlike the 

construction as described above, a construction may 
well be adopted which traces the two stages of a 
"critical probability calculation portion (goodness of 
fit quantification portion) that concretely calculates 

10 the value of the critical probability pi (x, y) , and a 

"weight calculation portion" that calculates the weight 
wl(p(x, y) ) on the basis of the output of the former 
portion. Anyway, the weight calculation portion 112 
calculates the weight wl(p(x, y) ) by using the variance 

15 a 2 estimated by the variance estimation portion 111, 

and the pixel values v(x) and v(y) . 

Besides, the pixel value calculation portion 113 
calculates the pixel value v' (x) by using the pixel 
values v(x) and v(y), and the weight wl(p(x, y) ) 

20 numerically calculated by the weight calculation 

portion 112. In other words, this calculation portion 
113 actually performs the processing of suppressing the 
noises of the original image, that is, the application 
of the coherent filter (hereinbelow, such an operation 

25 shall be expressed as "to apply the coherent filter") . 

In a case where, in the dynamic coherent filtering 
as stated above, the coherent filter is applied to a 



dynamic image constituted by K static images, the 
processing in the image processing unit 110 may well b 
such that all the static images are once reconstructed 
and are thereafter stored in the storage device 10M, 
and that the coherent filter is applied to the stored 
images as postprocessing later. However, the present 
invention is not restricted to such an aspect, but the 
processing of applying the coherent filter may well be 
executed in real time in the foregoing flow of the 
continuous scan, continuous projection data 
acquisition, continuous reconstruction and continuous 
display (hereinbelow, such processing shall be called 
"real-time coherent filtering") . 

In a preferred embodiment of the real-time 
coherent filtering, processing to be described below 
is executed each time a new image is obtained and 
reconstructed. On K static images having image Nos. M 
M - 1, •••/ and M - K + 1 , among a firstly obtained 
image (image No. 1) through the newest image (image 
No. M) , pixel values (scalar values) which pixels x of 
the same common points (same coordinates) have are 
arrayed to construct a K-dimensional vector value 
v(x) = (v M (x), v M _!(x), v M _ K+1 (x)). In this way, 

the coherent filter can be applied quite similarly to 
the "dynamic coherent filter processing" described 
before. Herein, the pixel value calculation portion 
113 does not actually compute all the elements of the 
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pixel value v' (x) , but it computes only the scalar 
value v^' (x) corresponding to the newest image (image 
No. M). As a result, a computing speed is enhanced, so 
that the newest image with noises suppressed can be 
5 displayed in real time. 

Another preferred embodiment of the "real-time 
coherent filtering" may well be so constructed that, 
when first K images have been obtained, the coherent 
filter is applied quite similarly to the above, so as 

10 to find V]/ (x) , • • •, and v K ' (x) , that a K-dimensional 

vector value is constructed of v(x) = (v M (x), 
v M-l' < x ) 9 v M-K+l' ( x ) ) fa y thereafter using K static 

images which have image Nos. M, M -1, •**, and 
M - K +1, and that the real-time coherent filtering 

15 is applied to the K-dimensional vector value. 

Incidentally, convenient is a construction in which, 
in the real-time coherent filtering, the dimensions K 
) of the pixel value vector v(x) are adapted to be 

occasionally altered by manual setting or automatic 

20 setting. 

(IV Case 2 of Applying Present invention to X-ray CT 
radiography; (To obtain Single image with Noises 
reduced) ) 

Next, there will be described a practicable 
25 example in which a single image with noises reduced is 

obtained unlike in the dynamic CT . Heretofore, 
generally in a case where a single CT image of little 
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noises is to be obtained, there have been considered 
for the accomplishment the two ways of a scheme wherein 
K CT images obtained by K times of iterative 
radiographing in which an ordinary X-ray tube current 
5 is used have their simple average taken in the 

direction of time (by way of example, the simple 
average is {v^ (x) + ••• + v K (x)}/K when the concept of 
the pixel value v(x) in this embodiment is utilized), 
with the result that the single image is obtained, and 

10 a scheme wherein the single image is obtained by one 

time of radiographing in which an X-ray tube current 
being K times as large as the ordinary one is used. In 
either of the schemes, the variance of noises becomes 
1/K times as large as in the single image which is 

15 obtained by one time of radiographing by the use of the 

ordinary X-ray tube current (for the reason that, in 
the former scheme, the noises are averaged, while in 
t the latter scheme, the quantity of X-irradiation 

and the variance of noises are in the inversely 

20 proportional relationship as stated before) . Hereto- 

fore, accordingly, the single CT image has been 
generally obtained by the latter scheme of one time of 
radiographing as long as the performance of an X-ray 
tube allows, in view of the efficiency of operations, 

25 or by the former scheme of K times of iterative 

radiographing in a case where the performance of an 
X-ray tube cannot withstand the large-quantity of X-ray 



tube current. 

With these schemes, however, the quantity of X- 
irradiation increases more as the noises are to be 
suppressed more. That is, the quantity of exposure of 
a patient P increases more. 

Here, when the coherent filter according to the 
present invention is utilized, the noises contained in 
such a single CT image can be suppressed without 
increasing the quantity of exposure of the patient P. 

First, in this embodiment, CT images numbering K 
are obtained by X-irradiation in a quantity smaller 
than an ordinary one, within a time period within which 
quite no change can be assumed as to a patient P. In 
this case, accordingly, as to a pixel x defined for the 
K CT images, a pixel value v(x) can be constructed as: 

v(x) = (v 1 (x), v 2 (x), v K (x)) 

• • • (3 mentioned before) 
as in the dynamic coherent filtering described before. 
By the way, in this case, the quantity of X-irradiation 
is small, and hence, the variance a 2 of noises 
contained in each individual CT image is comparatively 
large . 

Subsequently, a concrete form similar to the above 
equation (6) is given as a weight function in this 
case . 

Since, however, an image to be finally obtained is 
a single CT image in this case, it is considered to 
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employ a scheme as stated below. 

Employable is, for example, a scheme wherein, as 
in the dynamic CT radiography, the coherent filter is 
applied to all of the K images, and the resulting K 
5 images have their weighted average or simple average 

taken, so as to construct the desired single CT image, 
or a scheme wherein the coherent filter is applied to 
only a single image selected at will from among the K 
images (= the calculations in FIG. 4(c) are executed 

10 only for the single image selected) . 

Incidentally, according to the former scheme, a 
computation for obtaining the single averaged image 
from the K images which result from applying the 
coherent filter to all the K images is concretely 

15 given, for example, as follows: 

K 

v* f (x) = - £ v k'(x) • • • <H> 

K k = l 

Here, v*' (x) on the left side of Equation (11) denotes 
the pixel value (scalar) at the pixel x on the single 
20 image resulting from simple-averaging. On this 

occasion, the variance (a 2 )* of noises contained in 
v*' (x) , substantially becomes: 

-i2 

w(p(x, y) ) 



yeN(x) 



]T) w(p(x, y) ) 
yeN(x) 



(12) 



25 Besides, since w(p(x, y) ) > 0 holds here, a term 



summated on the right side of Equation (12) is always 
smaller than 1 (one) . Accordingly, a noise suppression 
effect based on this embodiment is superior to the 
noise suppression effect (a 2 /K) achieved by the 
prior-art scheme. Conversely speaking, if a noise 
suppression effect equivalent to that of the prior art 
is to be attained by employing this embodiment, the 
necessary quantity of X-irradiation may be smaller than 
in the prior-art scheme, and hence, the quantity of 
exposure of the patient P can be decreased. 

The effect attained by employing this embodiment 
is, for example, as shown in FIG. 7. FIG. 7(a) shows 
an image obtained by the prior-art scheme in which 
images based on K times of iterative radiographing are 
simply averaged in the direction of time, while 
FIG. 7(b) shows an image which is constructed by 
applying the coherent filter to all of K images and 
taking the simple average of the resulting K images, in 
accordance with this embodiment. It is understood that 
random noises seen in FIG. 7(a) are sufficiently 
suppressed in FIG. 7(b). (Incidentally, FIG. 7 shows 
the images obtained under the conditions that 
radiographing of 0.5-second scan was carried out three 
times (k = 3), that a rectangular area corresponding to 
13 x 13 pixels was used as the set N(x), and that the 
parameter C = 3 was set. 

Besides, FIG. 7(c) shows an image which is 



constructed considering the point that, since the noise 
suppression effect in FIG, 7(b) is remarkably 
excellent, the image becomes unfamiliar as compared 
with the image in the prior art, and which is 
concretely obtained by averaging FIGS. 7(a) and 7(b) . 
More specifically, letting the pixel value of the image 
in FIG. 7(c) be v**' (x) , and a pixel value in FIG. 7(a) 
be : 

v(x) (= {v x (x) + • + v K (x)}/K), 
the image in FIG. 7(c) is obtained by setting u = 1/2 
in : 

v**' (x) = u*v(x) + (1 - u)*v*' (x) (13) 
By the way, in the present invention, this embodiment 
itself (which can obtained the image in FIG. 7(b)) is 
not restrictive, but a construction is also preferable 
in which the value of u in the above equation (13) 
(subject to 0 < u < 1) , for example, is made adjustable 
according to the taste of a diagnostician. 
(V Case 3 of Applying Present invention to X-ray CT 
radiography; Multi-energy radiography (To reduce Noises 
of Plurality of sorts of images) ) 

Now, there will be described a case where the 
present invention is applied to so-called "multi-energy 
radiography". Here, the multi-energy radiography 
signifies a scheme wherein the energy distribution of 
irradiating X-rays is variously changed in such a way 
that a voltage which is applied to an X-ray tube 101 is 



altered, or that an X-ray filter which is constructed 
of a thin metal plate or the like is inserted between 
the X-ray tube 101 and a patient P, and wherein the 
identical patient P is radiographed under the plurality 
of energy distributions thus obtained. In this case, 
notwithstanding that the patient P does not change, 
images change depending upon the energy distributions 
of the X-rays, with the result that the plurality of 
sorts of images are obtained* 

Assuming here that the patient P has been 
radiographed K times while the energy distribution is 
being altered, a pixel value v(x) = (vi (x) , v 2 (x) , 
v K (x)) can be constructed by arraying scalar values in 
the pixel x of the K sorts of images obtained by the K 
times of radiographing. Accordingly, the coherent 
filter can be applied to the pixel value v(x) by 
utilizing the same means as in the dynamic coherent 
filtering in the section (III Case 1 of Applying 
Present invention to X-ray CT radiography; Dynamic CT 
radiography (To eliminate Noises from Plurality of 
images)), with the result that the K sorts of images 
with random noises suppressed can be obtained. 
(VI Cases of Applying Present invention to Medical 
image diagnostic equipment other than X-ray CT 
equipment) 

Now, there will be described embodiments in which 
the coherent filter according to the present invention 
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is applied to medical image diagnostic equipment other 
than X-ray CT equipment for example, nuclear magnetic 
resonance imaging equipment (so-called "MR I 
equipment" ) , nuclear medical diagnostic equipment 
("scintillation camera, SPECT equipment, or PET 
equipment"), and fluoroscopic equipment. 
(VI -1 Case of Applying Present invention to MR I 
equipment) 

In the first place, there will be described 
examples of applications • of the present invention 
concerning MR I equipment. Regarding the MR I equipment, 
the present invention is applicable to radiographic 
schemes hitherto called "averaging", "fast 
radiographing" and "three-dimensional radiographing" 
and also to a "radiographic scheme employing a 
plurality of pulse sequences", etc. 

First, the averaging is a radiographic scheme 
wherein, with the object of suppressing image noises by 
averaging as indicated by its name, radiographing is 
iteratively performed while no change appears in a 
patient P, and a plurality of (K) images obtained have 
their simple average taken, thereby to obtain a single 
image. When scalar values in the pixel x of the K 
images are arrayed, a pixel value v(x) = (vi (x) , 
V2 (x) , vj((x) ) can be constructed. Besides, 

regarding a weight function w, the above equation (5) 
may well be employed by way of example. Accordingly, 
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the coherent filter can be applied to the pixel value 
v(x) by utilizing quite the same means as in the 
dynamic coherent filtering in the section (III Case 1 
of Applying Present invention to X-ray CT radiography; 
5 Dynamic CT radiography (To eliminate Noises from 

Plurality of images)), with the result that a plurality 
of images with noises suppressed as are superior to 
the images based on the averaging can be obtained. 
Alternatively, a single image with noises suppressed 
10 can be obtained by utilizing quite the same processing 

as in the section (IV Case 2 of Applying Present 
invention to X-ray CT radiography; (To obtain Single 
image with Noises reduced) ) . 

Besides, the fast radiographing is a scheme 
15 wherein, when a change is appearing in a patient P, the 

patient P is iteratively radiographed, thereby to 
obtain a dynamic image which represents the change with 
) the passage of time, and it is also called "dynamic 

MRI". Also in this case, a plurality of (K) static 
20 images are similarly obtained, and hence, noises can be 

suppressed by applying the coherent filter in quite the 
same way as in the dynamic coherent filtering in the 
section (III Case 1 of Applying Present invention to 
X-ray CT radiography; Dynamic CT radiography (To 
25 eliminate Noises from Plurality of images) ) , as in 

the case of the averaging. Further, in the fast 
radiographing, it is also allowed to apply quite the 
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same processing as the real-time coherent filtering in 
the section (III Case 1 of Applying Present invention 
to X-ray CT radiography; Dynamic CT radiography (To 
eliminate Noises from Plurality of images)). 

The three-dimensional radiographing is a 
radiographic method which obtains a three-dimensional 
image at one time. Also in this radiographic method, 
the averaging (namely, the processing in which, with 
the object of suppressing noises, radiographing is 
iteratively performed so as to take the simple average 
of the resulting images) or the fast radiographing (in 
which, when a change is appearing in a patient P, the 
patient P is iteratively radiographed) is sometimes 
carried out. In such a case, processing is quite the 
same as in the above method of applying the coherent 
filter to the averaging or the fast radiographing, 
except the point that a pixel x is expressed as a 
three-dimensional vector x = (x, y, z) . Further, in 
order to apply the coherent filter when the three- 
dimensional radiographing is performed only once, a 
scheme which will be described in a section (VIII To 
Construct Pixel value v(x) from Single image) later is 
applicable . 

Meanwhile, in the MR I equipment, different 
physical parameters can be brought into images by 
altering a pulse sequence. The "pulse sequence" 
signifies the combination of the sequence of applying 
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RF pulses and a gradient magnetic field, intensities 
and time periods for applications, and a quiescent time 
period, and multifarious images can be obtained by 
changing the pulse sequence. As a very small number of 
5 examples of the pulse sequence which is typically used, 

there have been known ones which are called "T2*- 
weighted image", "Tl-weighted image" and "proton- 
density image" prefixed by the physical parameters that 
contribute most to the contrasts of the images. 

10 Anyway, notwithstanding that a patient undergoes no 

change, the different images are obtained depending 
upon the pulse sequences, with the result that a 
plurality of sorts of images are, in effect, obtained. 
Accordingly, the coherent filter can be applied in 

15 quite the same way as in the section (V Case 3 of 

Applying Present invention to X-ray CT radiography; To 
reduce Noises of Plurality of sorts of images) . 
) More concretely, assuming that the patient has 

been radiographed K times while the pulse sequence is 

20 being changed, a vector v(x) = (v x (x) , v 2 (x) , 

v K (x)) can be constructed by arraying scalar values in 
the pixel x of the images obtained by the K times of 
radiographing. Also in this case, accordingly, noises 
can be suppressed in such a way that the coherent 

25 filter is applied by utilizing quite the same means as 

in the dynamic coherent filtering, as in the case of 
the averaging. 
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(VI-2 Case of Applying Present invention to 
Scintillation camera, SPECT equipment or PET equipment) 

Secondly, there will be described an example of 
application of the present invention concerning a 
5 scintillation camera, SPECT equipment or PET equipment. 

Such an equipment turns the distribution of a 
pharmaceutic drug (radiopharmaceutic drug) in a subject 
into an image by giving the pharmaceutic drug, which 
contains a radioisotope (hereinbelow, termed M RI"), 

10 into the subject and then detecting gamma rays emitted 

by the natural decay of the RI . Especially in a case 
where the subject is a living body, the equipment is 
employed with the object of obtaining a so-called 
"functional map" which represents the degrees of the 

15 biochemical metabolism and blood circulation of the 

living body. The scintillation camera obtains the 
distribution the radiopharmaceutic drug as a two- 
) dimensional image as if the subject were examined by 

fluoroscopy. Besides, each of the SPECT equipment and 

20 the PET equipment obtains a three-dimensional image 

representative of the three-dimensional distribution of 
the radiopharmaceutic drug by utilizing computed 
tomography. 

In the scintillation camera, SPECT equipment or 
25 PET equipment, the quantity of gamma rays to be 

detected is usually very small, and hence, so-called 
"photon noise" ascribable to dispersion in the number 
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of photons as is caused by the quantum-mechanical 
property of the RI appears on an image. 

Heretofore, in order to suppress the photon noise 
without degrading the resolution of the image, there 
5 has been no other method than increasing the quantity 

of gamma rays to-be-detected. Concretely, there have 
been a method of increasing the quantity of the RI 
which is given, and a method of lengthening a 
radiographing time period. 

10 In a case where the quantity of the RI to be given 

can be increased (the increase of a dose is allowed 
for, for example, an experimental animal for a medical 
research) , the photon noise can be reduced by the 
increase, but such a method is attended with increase 

15 in the quantity of exposure of an operator. Besides, 

in a case where a human body is diagnosed (= the human 
body is the subject), the increase in the quantity of 
) the RI as stated above is not freely allowed originally 

(it is legally regulated) . Further, in case of 

20 radiographing with radiation which is emitted from the 

RI contained in an natural object, the quantity of the 
radiation cannot be increased. 

On the other hand, in a case where the method of 
lengthening the radiographing time period is to be 

25 employed, it is required for attaining a sufficient 

noise reduction effect that the subject is immobile for 
a very long time (for example, for several hours) . 
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This requirement is difficult of realization in a case 
where the living body (especially, the human body) is 
subject. Another problem has been that the 
processability (or cost performance) of the 
5 radiographing equipment lowers due to the long 

radiographing time period. In actuality, accordingly, 
it is only possible to set a suitable radiographing 
time period limited by an allowable time period (for 
example, within one hour) , and intense noise inevitably 

10 appears in the image. Therefore, the noise is 

suppressed only by the smoothing of the image. It has 
accordingly been considered unavoidable that the 
resolution of the image is spoilt as a drawback 
attendant upon the smoothing. 

15 In such a case, according to the present 

invention, processing as stated below can be executed. 
Within a time period equivalent to the suitable 
) radiographing time period mentioned above, a plurality 

of times of radiographing in shorter time periods are 

20 performed, thereby to obtain a plurality of images 

containing comparatively much noises. Subsequently, 
quite the same processing as in the section (IV Case 2 
of Applying Present invention to X-ray CT radiography; 
(To obtain Single image with Noises reduced) ) is 

25 applied, whereby an image with the noises suppressed 

more can be obtained without spoiling a spatial 
resolution. Moreover, the processing can be realized 
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without increasing the quantity of the RI . 

Further, also in the scintillation camera, SPECT 
equipment or PET equipment, there is a radiographic 
scheme in which, as in the X-ray CT or the like, 
5 an identical subject or patient is iteratively 

radiographed, and a change with the passage of time in 
the distribution of a radiopharmaceutic drug within the 
patient is traced, that is, a radiographic scheme which 
is called "dynamic scintigram", "dynamic SPECT" or 

10 "dynamic PET". Also in this case, quite the same 

processing as the dynamic coherent filtering or real- 
time coherent filtering in the section (III Case 1 of 
Applying Present invention to X-ray CT radiography; 
Dynamic CT radiography (To eliminate Noises from 

15 Plurality of images) ) is applicable, and the noises of 

images can be effectively suppressed. 

Further, in the PET equipment or SPECT equipment, 
there is sometimes used so-called "multiwindow 
radiography" wherein, by utilizing the fact that an RI 

20 emits gamma rays of energy peculiar to its sort, a 

plurality of sorts of RIs are simultaneously given to a 
subject or patient so as to obtain respective 
distribution images as to the plurality of sorts of RIs 
at one time. The multiwindow radiography is based on a 

25 technique wherein the energy of detected gamma rays is 

measured, and which of the plurality of sorts of RIs 
emitted the gamma rays is discriminated depending upon 
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if the levels of the energy fall within a numerical 
range (energy window) set artificially. 

More concretely, a suitable number of energy 
windows (for example, windows Wl and W2) are set, 
5 and images, such as an image based on gamma rays 

corresponding to one of the windows (for example, 
the window Wl) and an image based on gamma rays 
corresponding to the other window (for example, the 
window W2), are separately formed, whereby the 

10 plurality of RI distribution images stated above are 

obtained by one time of radiographing . 

That is, an identical subject is radiographed 
using an imaging equipment, and the plurality of sorts 
of images are obtained under different treatment 

15 conditions. More generally speaking, the multiwindow 

radiography is a "radiographic scheme wherein, using K 
sorts of RIs, K sorts of images different from one 
another are obtained by one time of radiographing". 
Accordingly, the coherent filter can be applied in 

20 quite the same way as in the section (V Case 3 of 

Applying Present invention to X-ray CT radiography; To 
reduce Noises of Plurality of sorts of images) . 

Thus, also in this case, as to a pixel x on such K 
sorts of images, its pixel value v(x) = (v^ (x) , 

25 V2 (x) , v K (x)) can be constructed. Besides, 

although the variance of noises contained in the pixel 
value v(x) differs depending upon the position of the 
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pixel x and the sort of the image, it may be considered 
substantially proportional to the total amount of gamma 
rays detected every sort of the image, and it can be 
estimated from this property. As thus far described, 
the coherent filter can be applied also to the 
"multiwindow radiography" by applying the present 
invention, with the result that K sorts of images with 
random noises suppressed can be obtained. 

By the way, in a nuclear medical diagnostic 
equipment, there are reconstruction methods effective 
for enhancing a resolution (an MU-EN reconstruction 
method for a fan beam, a superresolution reconstruction 
method for a fan-beam or parallel-hole collimator, 
etc.). It is not restricted in principle at all to 
combine these methods and the coherent filter according 
to the present invention. It is consequently possible 
to create images which have high resolutions, and in 
which random noises are sufficiently suppressed. 
(VI-3 Case of Applying Present invention to 
Fluoroscopic equipment) 

Thirdly, there will be described an example of 
application of the present invention concerning a 
fluoroscopic equipment. Supposed as the "fluoroscopic 
equipment" termed here is an equipment wherein, while a 
patient is being continuously exposed to X-rays of low 
dose, it is radiographed iteratively every predeter- 
mined time period (for example, 1/24 second), and 
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wherein X-ray images obtained for individual "frames" 
each of which indicates one time of radiographing 
are successively displayed, whereby a so-called 
"fluoroscopic image (dynamic image) " can be observed. 
In general, accordingly, an equipment which includes a 
construction called M C arm" or arm" corresponds to 
the fluoroscopic equipment. Besides, an aspect to be 
stated below is naturally applicable to the foregoing 
X-ray CT equipment 100 as long as the functions stated 
above are installed. 

In such a fluoroscopic equipment, the fluoroscopic 
image contains much random noises ascribable to photon 
nose because of the low dose mentioned above. When the 
quantity of X-irradiation is sufficiently increased in 
order to solve such a problem, the quantity of exposure 
increases drastically, not only in the patient, but 
also in an inspector, for the reason that the 
continuous X-irradiation is performed. 

When the coherent filter according to the present 
invention is applied, it is permitted to solve these 
problems and to offer a fluoroscopic image with the 
noises suppressed satisfactorily. 

Concretely, image Nos. 1, 2, and M shall be 

successively assigned to the images obtained for the 
respective frames, for the sake of the ensuing 
description. The newest K images, namely, the images 
of image Nos. M, M - 1, and M - K + 1 are stored 
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in a memory, and these K images are subjected to quite 
the same processing as the dynamic coherent filtering 
or real-time coherent filtering in the section (III 
Case 1 of Applying Present invention to X-ray CT 
5 radiography; Dynamic CT radiography (To eliminate 

Noises from Plurality of images) ) . 

Incidentally, apart from the medical image 
diagnostic equipment, the above method is applicable 
also to "fluoroscopic equipment" used for so-called 

10 "nondestructive testing". 

Besides, the present invention is similarly 
applicable to a tomogram which is obtained by an 
ultrasonic diagnostic equipment. Also in this case, it 
holds true that images are successively generated while 

15 a subject is being continuously radiographed. It is 

accordingly possible to apply the coherent filter by 
the same construction as in the fluoroscopic equipment 
% and to suppress noises. 

(VII Case of Applying Present invention More 

20 generally; Color image, etc.) 

The present invention is also applicable to 
"general equipment" handling images, other than the 
medical image diagnostic equipment stated above, or 
"general images". There will now be explained examples 

25 of applications of the present invention to such 

general equipment or images. Incidentally, a case of 
handling general "color images" as is an especially 
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important field shall be chiefly described here. 

In a color image, a three-dimensional vector whose 
elements are the brightnesses of the three primary 
colors of light can be constructed as to each pixel x 
constituting the image, as already stated. That is: 

v(x) = (v R (x), v G (x), v B (x)) (3. 1) 

(where vr, vq and v B signify the lights of R (red) , G 
(green) and B (blue) being the three primary colors of 
light) . As seen from FIG. 2, this corresponds to the 
case of "constructing pixel values v(x) from a 
plurality of sorts of images''. 

Besides, in a case where the color image is taken 
by, for example, a digital camera which is constructed 
of a CCD imaging device, a CMOS imaging device or the 
like, the variance of noises contained in the pixel 
values of each pixel is determined by the quantity of 
light having entered the pixel x. Concretely, the 
scalar values v R (x), v G (x) and v B (x) and the variances 
a 2 R (x), a 2 G (x) and a 2 3(x) of noises contained in the 
scalar values are in substantially proportional 
relationships, respectively. The estimative values of 
the variances of the noises can be found by utilizing 
this property. Accordingly, the concrete form of the 
weight function w as in the above equation (6) is given 
and is used for constructing the coherent filter of the 
above equation (4), whereby pixel values v' (x) = 
(v' R (x), v' G (x), v'b(x)) with the noises suppressed can 
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be obtained. 

When the word "color" is grasped in the broadest 
sense so as to comprehend, for example, infrared rays 
and ultraviolet rays in addition to the red, green and 
blue lights mentioned above, the equation (3. 1) can 
also be constructed as a vector of higher order or as a 
vector independent of the colors R, G and B. By way of 
example, it is known that the method of representing a 
color image by using the pixel values of vector values 
is not always limited to the three primary colors only, 
and the present invention may well be applied to a case 
of employing such a vector representation not using the 
three primary colors. Besides, cases of four primary 
colors including ultraviolet light Vjj, five primary 
colors with violet Vp separated (incidentally, this 
corresponds to the vision of butterflies), etc. are 
considered. In the case of, for example, the five 
primary colors, the pixel x is handled as having the 
five-dimensional vector values of V(x) = (Vr(x), V G (x), 
V B (x), Vu(x), Vp(x)), whereby the coherent filter can 
be applied as in the foregoing construction. 

Further, in a photographing equipment wherein, in 
order to measure, for example, the light reflection 
spectrum of a subject, the subject is photographed by a 
monochromatic camera iteratively at individual 
illumination wavelengths while the wavelength of 
illumination light for illuminating the subject is 
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being changed-over stepwise or continuously, or the 
subject is photographed iteratively by changing-over 
color filters stepwise, thereby to perform the 
photographing at a multiplicity of wavelengths, the 
5 pixel values v(x) can be constructed on the basis of a 

large number of images (each being of single wavelength 
or single color) . That is, the coherent filter can be 
applied by utilizing quite the same method as in the 
section (V Case 3 of Applying Present invention to X- 
10 ray CT radiography; To reduce Noises of Plurality of 

sorts of images), with the result that a plurality of 
sorts of images with random noises reduced are 
obtained. 

Still further, many of such imaging equipment 
15 usually carry out one time of photographing in order to 

obtain a single image. However, when the present 
invention is applied, a plurality of (K) images should 
better be obtained by iterating photographing of short 
time period within the same photographing time period, 
20 as stated in the section (IV Case 2 of Applying 

Present invention to X-ray CT radiography; (To obtain 
Single image with Noises reduced) ) concerning the "X- 
ray CT radiography". The reason therefor is that 
(although the K images contain much random noises 
25 because the quantity of light becomes small due to the 

short-time photographing, ) a single image with the 
random noises suppressed more intensely can be 
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constructed by applying the coherent filter to the K 
images. Moreover, in a case where the K images are, 
for example, color images, they may well be handled as 
images "in a plurality of combinations (K combinations) 
5 each of which consists of a plurality of sorts (three 

sorts) of images obtained in correspondence with the 
three primary colors of light''. More specifically, 
a three-dimensional vector whose elements are the 
brightnesses of the three primary colors of light can 

10 be constructed every pixel x constituting each image, 

so that the pixel values v(x) of a 3K-dimensional 
vector can be constructed for all of the K images. The 
coherent filter may well be applied to the pixel values 
of the 3K-dimensional vector. 

15 Besides, in an imaging equipment, such as video 

camera, which takes a dynamic image, it is possible to 
apply the dynamic coherent filtering or real-time 
% coherent filtering stated in the section (III Case 1 of 

Applying Present invention to X-ray CT radiography; 

2 0 Dynamic CT radiography (To eliminate Noises from 

Plurality of images) ) . Further, the coherent filter 
can be applied by handling the dynamic image as images 
"in a plurality of combinations (K combinations) each 
of which consists of a plurality of sorts (three sorts) 

25 of images obtained in correspondence with the three 

primary colors of light" as in the foregoing. 

Incidentally, regarding images taken by a video 
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camera, a noctovision camera, a fast photographing 
camera, etc. (generally, termed "imaging equipment" 
below) other than the digital camera stated above, the 
present invention is applicable in the same manner as 
5 described above. Especially, the noctovision camera or 

the fast photographing camera is a camera which is 
originally used in a situation where the quantity of 
light per frame is small, and images contain much 
random noises. It is unfavorable, however, to spoil a 

10 spatial resolution, and to spoil a temporal resolution 

particularly in the fast photographing camera. 
Therefore, the noctovision camera or the fast 
photographing camera is one of equipment most suited to 
apply the present invention. 

15 Also in the imaging equipment, accordingly, images 

with noises suppressed more can be obtained by such 
photographing method. 

Meanwhile, in an equipment, for example, video 
tape recorder or DVD player, which does not execute 

2 0 w imaging", but which is capable of reproducing an image 

that is based on a signal recorded on any recording 
medium or transmitted, when the coherent filter 
according to the present invention is applied to the 
image to-be-reproduced, a reproduced image with random 

25 noises suppressed can be obtained. 

Further, in a case where a static image (so-called 
"pause image" or image for printing) is generated and 
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displayed in such an equipment, a static image with 
random noises suppressed is obtained by applying the 
coherent filter, and it is also suitable for the 
printing, etc. For this purpose, it is possible to 
5 employ the method stated in the section (III Case 1 of 

Applying Present invention to X-ray CT radiography; 
Dynamic CT radiography (To eliminate Noises from 
Plurality of images) ) , the method stated in the section 
(V Case 3 of Applying Present invention to X-ray CT 
10 radiography; To reduce Noises of Plurality of sorts of 

images) , a method to be stated below in a section (VIII 
To Construct Pixel value v(x) from Single image), or 
the combination of these methods. 

(VIII To Construct Pixel value v(x) from Single image) 
15 The examples of applications thus far described 

correspond to those examples in the classification 
shown in FIG. 2 in which the pixel values v(x) are 
respectively constructed "from a plurality of images'' 
and "from a plurality of sorts of images". In the 
20 present invention, in addition to the two methods of 

constructing the pixel values v(x), to "construct a 
pixel value v(x) from a single image" is considered as 
also shown in FIG. 2, and this method of construction 
can be performed as stated below. 
25 In a single two-dimensional image, for example, a 

single CT image acquired by the X-ray CT equipment 100, 
the MRI equipment or the like, or a single image based 



on the monochromatic photographing by the imaging 
equipment or the like, the value of a pixel x is 
constructed of only one scalar value. In other words, 
a pixel value v(x) being a vector value becomes a 
unidimensional vector . 

One method for applying the coherent filter 
according to the present invention, in such a case, is 
to consider the pixel value v(x) as the "unidimensional 
vector'' itself. That is, in the above equation (3), 
K = 1 is put to obtain: 

v(x) = v(x) • • • (3. 2) 

Thenceforth, the construction of the coherent filter or 
the suppression of image noises can be executed in the 
same way as in the foregoing procedure. The present 
invention covers also such a case where the "pixel 
value" v(x) being a "vector value" is the 
unidimensional vector, namely, the scalar value. 

Besides, a weight w3 in this case can be 
calculated using, for example, the following equation: 

w3 (p (x, y) ) = 



Here, D has a property similar to that of the 
parameter C in Equation (6). That is, when D is small, 
the weight w3(p(x, y) ) becomes smaller, and when D is 
large, the weight becomes larger. 




: for | v(x) - v(y) | < Da 



••• (14) 



0 



: for | v(x) - v(y) | > Da 
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Images which are based on pixel values v' (x) which 
are calculated from such a weight w3(p(x, y) ) and the 
above equation (4), become as shown in FIG. 8 by way of 
example. FIGS. 8(a), (b) , (c) and (d) are enhanced 
5 images each showing an example of 3D sampling of blood 

vessels, and FIG. 8 (e) is an image obtained by sampling 
the gray matter of a brain as another sampling example. 
FIG. 8(a) shows a CT image (original image) which is to 
be processed by the coherent filter. FIGS. 8(b) and 

10 (d) show results obtained in such a way that the 

coherent filter is applied using a rectangular area of 
3x3 pixels as the set N(x) in Equation (6) . The 
result in FIG. 8(b) corresponds to a case where D = 9 
is set in the above equation (14) , while the result in 

15 FIG. 8(d) is such that the image in (c) obtained by 

applying the coherent filter has been further subjected 
to threshold processing. Besides, shown in FIG. 8(c) 
is the image which is such that the original image in 
(a) has been subjected to threshold processing. 

20 When the images of FIGS. 8(a) and 8(b) are 

compared, it is understood that, owing to the 
application of the coherent filter, the image of 
FIG. 8(b) has its noises properly suppressed in 
comparison with that of FIG. 8(a), so a quality 

25 satisfactory for a diagnosis is guaranteed. 

When the images of FIGS. 8(c) and 8(d) are 
compared, it is understood that the image of FIG. 8(c) 
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obtained by subjecting the original image to only the 
threshold processing is difficult of blood vessel 
sampling based on a threshold value, whereas the image 
of FIG. 8(d) obtained by the threshold processing after 
5 the application of the coherent filter permits the 

blood vessel sampling based on the threshold value, 
without omitting even fine blood vessels. 

Besides, according to FIG. 8(e), it is understood 
that the region of the gray matter has been sampled. 

10 According to such a scheme, the effects as shown 

in FIG. 8 are really attained. Since, however, the 
unidimensional vector value v(x) = v(x) is used for 
measuring the goodness of fit between the pixels x and 
y, the performance of discrimination for the goodness 

15 of fit is comparatively low, in other words, only the 

critical probability p(x, y) which is statistically 
uncertain is obtained. Therefore, the effect of 
suppressing random noises without spoiling a resolution 
cannot be said considerably high. In such a case, 

20 accordingly, another scheme to be stated below should 

be adopted more preferably. Effective is the scheme 
wherein a multidimensional vector is constructed using 
a set Z(x) which consists of some pixels in the 
vicinity of the pixel x, so as to employ the 

25 multidimensional vector as the pixel value v(x) of the 

pixel x. By way of example, as shown in FIG. 9, let's 
consider a five-dimensional vector wherein, using the 
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set Z(x) which consists of four pixels (x + 1, y), 
(x - 1, y) , (x, y + 1) and (x, y - 1) adjoining the 
pixel x = (x, y) around it (hereinbelow, the pixels 
which are the elements of the set Z (x) around the 
pixel x shall be expressed as "z 1 (x) , z 2 (x) , z 3 (x) and 
Z4(x)"), the scalar values of the four pixels are added 
to the scalar value of the pixel x itself, namely, 
v(x) = v(x, y) . That is: 

v(x) = (v(x, y) , v(x +1, y) , v(x - 1, y) , 
v(x, y + 1), v(x, y - 1) (3. 3) 

Such vectors can be similarly constructed as to all 
pixels which constitute an image G shown in FIG. 9 
(incidentally, FIG. 9 shows only a part of the whole 
image on an enlarged scale) . 

In this way, the coherent filter according to the 
present invention is constructed using the pixel value 
v(x) of the multidimensional vector constructed from 
the single image, whereby the goodness of fit can be 
decided at a higher precision than in the case of 
employing the pixel value of the unidimensional vector. 
As a result, the random noises can be suppressed more 
intensely without spoiling the resolution of the image. 

By the way, in the case where the pixel value is 
thus constructed, the foregoing null hypothesis H for 
use in quantifying the goodness of fit, "the pixels x 
and y have the same pixel values when their noises are 
eliminated 7 ' signifies nothing but the null hypothesis 
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H' that "the pixel x and the surrounding pixels z 1 (x) - 
z 4 (x), totaling five, and the pixel y and the 
surrounding pixels z± (y) - z 4 (y) , totaling five, have 
the same pixel values as those of the corresponding 
5 pixels except the noises, respectively". 

More concretely, in a case where, as shown in 
FIG. 9, a certain pixel x and a certain pixel y which 
exists within a set Z (x) (this is, for example, an area 
around the pixel x as stated above) are both existent 

10 inside the region El of any object El which appears on 

an image G as shown in the figure by way of example, 
the value of the critical probability p(x, y) computed 
by, for example, Equation (9) becomes small, and the 
value of the weight w(p(x, y) ) computed by, for 

15 example, Equation (6) becomes large. The reason 

therefor is that, in this case, pixels z 2 (x) and z 4 (x) 
and pixels z 2 (y) and z 4 (y) lie in the region of the 
image El of the particular object, while pixels z 1 (x) 
and z 3 (x) and pixels z 1 (y) and z 3 (y) lie in the image 

20 of an object E2, so vector values V(x) and V(y) 

expressed by Equation (3. 3) are similar. 

Now, in the above embodiment, when note is 
especially taken of the pixel x which lies at the edge 
of the image El of the object, it is conjectured that 

25 pixels y as to which the critical probability p(x, y) 

in the case of abandoning the null hypothesis H' is low 
will not usually be found in a very large number (in 



FIG. 9 by way of example, regarding even the pixel x 
and a pixel y p both of which exist inside the edge E1E 
of the image of the object El, the critical probability 
p(x, y) in the case of abandoning the null hypothesis 
H' becomes higher as compared with the critical 
probability concerning the pixels x and y shown in the 
figure) . 

More generally speaking, regarding a pixel x 
existent at the edge of the image of a certain object, 
and most of yeN(x), the value of the critical 
probability p(x, y) computed by, for example, Equation 
(9) becomes large, and the weight w(p(x, y) ) computed 
by, for example, Equation (6) becomes small. Then, a 
smoothing effect is feeble even by weighted averaging, 
so. that the noise suppression effect based on the 
coherent filter as computed by, for example, 
Equation (4) is scanty. 

Therefore, another scheme with the above scheme 
further improved will be described below. In this 
scheme, a pixel value v(x) is constructed as: 

v(x) = (v(x, y) , v x , v 2 , v 3 , v 4 ) (3. 4) 

where v^ - V4 are such that v(x +1, y) , v(x - 1, y) , 
v(x, y + 1) and v(x, y - 1) in the above equation 
(3. 3) are rearrayed in conformity with a predetermined 
rule (for example, in the order of larger values) . 

In this way, the null hypothesis H' turns into the 
more moderate null hypothesis H' ' that "the pixel x and 



the surrounding pixels (x) - Z4 (x) , totaling five, 
and the corresponding pixel y and the surrounding 
pixels zi (y) - Z4 (y) , totaling five, have the same 
pixel values except noises [while allowing the changes 
of sequences within the pixels z^ (x) - Z4 (x) and within 
the pixels z^ (y) - Z4 (y) ] , respectively (the contents 
in the brackets [ ] have been added) . 

Besides, in this case, in the relationship between 
pixels x and y which exist inside the edge E2E of the 
image E2 of an object as shown in FIG. 10 by way of 
example, the critical probability p(x, y) in the case 
of abandoning the null hypothesis H' ' becomes low 
(because the pixels Z£ (x) , Z3 (x) and Z4 (x) and the 
pixels zi(y), Z2 (y) and Z4 (y) agree as v^ through V3 in 
the above equation (3. 4), and the pixel z^ (x) and the 
pixel Z3 (y) similarly agree as V4) , with the result 
that a noise suppression effect becomes greater than in 
the scheme based on the construction of the pixel value 
v(x) conforming to the above equation (3. 3), and a 
spatial resolution is not spoilt. 

As described above, in the present invention, 
a pixel value v(x) being a vector value can be 
constructed from a single image, so that noises can be 
suppressed without spoiling a spatial resolution, by 
employing the coherent filter. Since such a scheme 
suffices with the single image, it is generally 
applicable to any images fundamentally. That is, this 
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scheme is applicable, not only to the foregoing various 
examples of applications (X-ray CT equipment, medical 
image diagnostic equipment, color images, etc.), but 
also to all images. Besides, even in a case where the 
5 image is a dynamic image (which can be seen as a 

plurality of static images as stated above several 
times), this scheme (in which "the pixel value v(x) is 
constructed from the single image") may, of course, be 
applied to each of the static images. Further, in case 

10 of employing this scheme for a three-dimensional or 

multidimensional image, a set Z (x) can be constructed 
of, for example, a pixel Xl itself and six pixels lying 
in contact with the pixel Xl • Such expanding methods 
can be suggested very easily, and are self-explanatory. 

15 By the way, in the present invention, the coherent 

filter may, of course, be constructed (as already 
stated) by using the pixel value v(x) which is 
\ constructed in accordance with each of the individual 

schemes "from a plurality of images", "from a plurality 

20 of sorts of images" and "from a single image" as shown 

in FIG. 2, but the coherent filter can also be 
constructed on the basis of a pixel value v(x) which is 
constructed by properly combining the respective 
schemes and which is expressed as a vector value of 

25 larger dimensions. 

As already stated in the section VII, assuming by 
way of example that there are two color images 
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concerning an identical subject, as to a certain pixel 

scalar values v R (x), vq(x) and V3(x) indicated on 
the right side of the above equation (2. 2) are derived 
from the two images as the pixel value v(x) of the 
5 pixel x as is a six-dimensional vector: 

v(x) = (v R1 (x) , v G1 (x) , v B1 (x), v R2 (x), v G2 (x), 
v B2 (x)) (3. 5) 

When the pixel value v(x) expressed by a vector of 
larger dimensions is used in this manner, a more 

10 reliable value can be calculated as the critical 

probability p(x, y) which is computed by, for example, 
Equations (9) and (6) . It is consequently permitted to 
suppress random noises more intensely without spoiling 
the resolution of the image. 

15 (IX Most general form of Coherent filter in 

Embodiment) 

Now, apart from the above description along 
FIG. 2, a more general form concerning the coherent 
filter introduced as in the above equation (4) will be 

20 described. 

In the general form (Equation (4) ) of the coherent 
filter, or the various examples of applications, there 
has been employed the null hypothesis H that "the 
pixels x and y have the same values when their noises 

25 have been removed" or that ^v(x) = v(y) except 

differences ascribable to the noises of both the 
pixels" stated at the beginning, but this null 
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hypothesis can be given a more general expression as 
stated below. 

A model function f (a, v(y)) which includes an 
unknown M-dimensional parameter a = (a^, a^) 
5 (where M < K) is introduced, whereupon a null 

hypothesis Hq "v(x) = f (a, v(y)) + £ (where £ conforms 
to a (known) probability distribution)" is built up as 
to two certain vectors x and y which represent the 
positions of pixels. On this occasion, the estimative 

10 value a~ of the parameter a is determined from v(x) and 

v(y) by employing a suitable fitting method (for 
example, the least squares method to be explained 
later) , whereby the above null hypothesis can be said 
in other words as its equivalent proposition Hq "v(x) = 

15 f(a~, v(y)) + £, (where £ conforms to a (known) 

probability distribution)". Here, f (a~, v(y)) 
mentioned above is such that the pixel value v(y) of 
the pixel y is transformed so as to have the highest 
goodness of fit with the pixel value v(x) of the pixel 

20 x by optimally adjusting the degree of freedom (namely, 

the parameter a) allowed by a model which the function 
f describes. Besides, such an optimum parameter is a". 

As a result, a critical probability p(x, y) is 
permitted to be concretely computed using the known 

25 probability distribution (with which £ is supposed to 

conform if the null hypothesis is correct), so that a 
weight w(p(x, y) ) can be computed. Therefore, a 
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weighted average v') c (x) can be computed by: 



£ w(p(x, y) )f (a~, v(y) ) 




(15) 



10 



15 

\ 

20 



This is the most general form of the coherent filter in 
this embodiment. Incidentally, when the above equation 
(4) and this equation (15) are compared, it is 
understood that the former is nothing but a more 
concrete form obtained by actually putting the 
following in the latter: 



(IX-1 Example of application based on Above equation 
(15) ; Construction of Time-density curve based on 
Dynamic image) 

There will be stated another example of 
application of the present invention employing the form 
of the coherent filter as indicated by Equation (15) . 
Here, the example of application is an embodiment in 
which, in quantitatively measuring a time-density curve 
that handles the dynamic CT image obtained by the 
medical image diagnostic equipment, the time-density 
curve with random noises sufficiently suppressed is 
intended to be acquired by applying the present 
invention . 

First, the "dynamic image" termed here is the 
combinations between a series of images (image Nos. 1, 
2, K) obtained by iteratively radiographing an 

identical subject, and the radiographing times 



f (a^, v (y) ) = v (y) 



(16) 
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(ti, tp() of the respective images. It is, for 

example, the dynamic CT image obtained by the medical 
image diagnostic equipment, namely, the X-ray CT 
equipment 100, or it is radiographed by employing the 
5 fluoroscopic equipment, scintillation camera, MR I 

equipment, SPECT equipment or PET equipment, ultrasonic 
diagnostic equipment, or the like. In addition, the 
dynamic image may well be composed of a series of 
images obtained by iterative photographing with the 

10 general photographing equipment, namely, the digital 

camera, video camera, noctovision camera, fast 
photographing camera, or the like. 

Besides, the "time-density curve" termed here is a 
curve which represents a change with the passage of 

15 time in the image density value of a specified part 

within the dynamic image. Especially in the medical 
image diagnostic equipment, for the purpose of 
examining the details of, e.g., a blood flow movement, 
a metabolic function in a human body tissue or the 

20 like, the change with the passage of time in, e.g., a 

contrast medium density in the specified tissue of the 
human body is measured as the time-density curve. 
Besides, in an astronomical observation or the like, 
for the purpose of analyzing, e.g., a change in the 

25 brightness of a specified celestial body, the time- 

density curve is utilized. When more formally stated, 
the "time-density curve" is expressed as a sequence of 
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pairs {<t) c , dfc> (k = 1, 2, K) } where denotes 

the density value of a certain part at a time t^. 
Besides, in many applications of the time-density 
curve, the absolute value of the density value d^ is 
5 not always necessary, but it suffices to obtain only an 

increment (d^ - d^) taken with reference to the first 
image 1. Further, in many of such applications, it 
suffices to obtain only data A(dk - d^) proportional to 
the increment (d^ - d^) (where A denotes an unknown 
10 proportionality coefficient) . In this case, 

accordingly, a sequence of pairs {<tj c/ Afd^ - d^)> (k = 
1, 2, ••*, K) } is the time-density curve to-be- 
acquired. 

In order to acquire such a time-density curve, the 
15 sequence of pairs {<t} c/ vj c (x)>} or {<tj c , A(vj c (x) - 

vi(x)>} may, in principle, be constructed using the 
scalar value v^fx) of a pixel x included in the part as 
to which the time-density curve is to be measured, in 
each of the images k (k = 1, 2, K) constituting 

20 the dynamic image. 

In practical use, however, there is the problem 
that, since random noises are contained in the dynamic 
image obtained by the medical image diagnostic 
equipment or the like, the time-density curve to be 
25 originally measured cannot be precisely acquired. 

Further, in the practical use, a so-called 
"partial volume effect" occurs in the dynamic image. 



The partial volume effect is the phenomenon that the 
image of a minute object within the subject or patient 
is expressed by a small number of pixels on the dynamic 
image, and that the small number of pixels are 
influenced by the image of adjacent objects in the 
subject, so the pixel values of the small number of 
pixels exhibit only small fluctuations (though they are 
proportional to the fluctuations of density values 
which are to be originally measured) . In other words, 
the pixel values of the small number of pixels contain 
only slight signals. Accordingly, in the case where 
the partial volume effect occurs, the sequence of pairs 
{<tj c , v) c (x)>} (k = 1, 2, • * • , K) has a very low signal 
level at any pixel x, and it is influenced by the 
change of a density value in a tissue which is not to 
be originally measured. Further, random noises are 
existent. These result in the problem that the time- 
density curve {<tk, d] c >} to be originally measured 
cannot be precisely acquired. 

In order to suppress the random noises, therefore, 
a method based on temporal averaging has hitherto been 
employed as stated below. First, the sequence of pairs 
{<tj c , vj c (x)> (k = 1, 2, K) } is constructed 

concerning the pixel x which is included in the part 
to-be-measured. Subsequently, a temporal average value 
vj c #(x) is computed by, for example, the following: 



80 



n 



Z v k + j(x)Gj 



v k #(x) = 3 m (17) 

Z^j 
j=m 

Besides, the temporal average value of times 

corresponding to the above temporal average value is 

computed by the following: 

n 

Z tk + jMGj 

5 t k #(x) = j = m ••• (18) 

Z^j 

j = m 

(where m < n denotes suitable integers, and Gj denotes 
a suitable weighting factor) . Thereafter, a time- 
density curve {<tj c #, v k #(x)> (k = 1, 2, K) } is 
constructed and employed. 

10 Alternatively, in order to suppress the random 

noises, a method based on spatial averaging has 
hitherto been employed as stated below. A set R of 
pixels as substantially corresponds to the part to-be- 
measured (that is, an ROI (Region of Interest) on the 

15 dynamic image) is set. In the image of No. k, the 

spatial average value V] C (R) of the scalar values v^fx) 
of all pixels xgR included in the set is calculated by 
Equation (19) given below. Thereafter, the sequence of 
pairs {<t] c , v^(R)> (k = 1, 2, K) } is constructed 

20 and is employed as the time-density curve. 
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v k (R) = - ■ • • (19) 

xgR 

Further, there has been employed a scheme wherein 
the random noises are more suppressed by utilizing both 
the temporal averaging and the spatial averaging. 
5 These prior-art schemes, however, have had the 

problems that, when the temporal averaging is executed, 
a temporal resolution is spoilt, and that, when the 
spatial averaging is executed, the changes with the 
passage of time in the densities of parts other than 

10 the part to be originally measured mix into a 

measurement value . 

In order to solve such problems and to acquire a 
preciser time-density curve, the coherent filter 
according to the present invention can be applied on 

15 the basis of the section (IX Most general form of 

Coherent filter in Embodiment) . 

First, there will be explained a null hypothesis 
which is to be utilized in the coherent filter of this 
embodiment. Assuming that the true time-density curve 

20 at the part to be measured is {<t) c , d^> (k = 1, 2, 

K)}, in a case where {<t] c , A(d k - d^) > (k = 1, 2, 
K) } (where A denotes an unknown coefficient) being the 
linear transformation of the true time-density curve is 
intended to be measured, a set R of pixels as substan- 

25 tially corresponds to a part to-be-measured is set* 

If, as to any pixel xeR being the element of the set R, 
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a condition Q: "the pixel x reflects the true time- 
density curve favorably and is hardly influenced by 
the density changes with the passage of time in the 
other parts" is fulfilled, it can be assumed that 
Equation (20) given below holds by considering the 
influences of the partial volume effect and random 
noises as to a pixel value (being a vector value) 
v(x) = (v 1 (x) / v 2 (x), v K (x)): 
v k (x) = p(x) d k + q(x) + y k (x) 

(k = 1, 2, • • K) • • • (20) 

Here, p(x) and q(x) denote unknown coefficients which 
differ every pixel x, but which do not change depending 
upon the image No. k (that is, the radiographing time 
t k ) , and they are obtained by modeling the partial 
volume effect. Besides, Y k (x) is obtained by modeling 
the random noises. Although the value of Y k (x) differs 
every pixel x and every image No. k, the expected value 
thereof is 0 (zero), and the statistic distribution 
thereof depends upon neither the pixel x nor the image 
No. k. (By the way, in the ensuing description, a 
case where the statistic distribution is a normal 
distribution of average 0 and variance a 2 shall be 
exemplified for the sake of brevity, but it is only 
required that the statistic distribution is substan- 
tially approximated by any known distribution, and how 
to modify the following embodiment so as to be suited 
to that case is obvious.) 
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Under the above assumptions, if the proposition 
that "both the pixels x and y satisfy the condition Q 
(mentioned above)" holds regarding any two pixels x and 
y which are the elements of the set R, it can be 
5 demonstrated that the relationship of the following 

equation holds: 

v k (x) = aiVfcty) + a 2 + £k 

(k = 1, 2, • • K) . • • (21) 

Here, a^ and a 2 denote unknown coefficients which 
10 differ every combination x, y of the pixels, but which 

do not change depending upon the image No. k (that is, 
the radiographing time t^) . Besides, denotes the 
random noises, and the value thereof differs every 
pixel combination x, y and every image No., but the 
15 expected value thereof is 0. 

Equation (21) is derived as stated below. When 
an equation obtained by substituting y into x in 
Equation (19) : 

V]c(Y) = P(y) d k + q(y) + yy^(y) (22) 
20 is reduced, the following is obtained: 

= v k (y) - q(y) - y k (y) ... 

p (y) 

By substituting this equation into Equation (20) , the 
following is obtained: 



25 



v k (x) = ^ v k (y) + | q(x) - ^ q(y) 

p(y) I P(y) 



••• (23) 

Accordingly, Equation (21) is derived by putting: 
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P(x) 
ai = 

p(y) 

a 2 = q(x) - 



(24) 



P(Y) 

£k = Yk< x ) ~ Yk(Y) 
P(Y) 

Here, and a 2 in Equation (24) denote parameters 

which represent the partial volume effect, and £ k in 
the same denotes random noises. 
5 It has been indicated from the above that the 

proposition that "the pixels x and y both satisfy the 
condition Q" is equivalent to the null hypothesis Hq' 
that "v k (x) = a x v k (y) + a 2 + £ k (k = 1, K) holds" 

Next, there will be explained a method for 

10 transforming the null hypothesis Hq' that "v k (x) = 

a l v k(y) + a 2 + £k (k = 1, K) holds", into a 

proposition in a form which is, in effect, equivalent 
thereto and which can be actually tested. When this 
null hypothesis is stated in a mathematically strict 
% 15 expression anew, it becomes the null hypothesis Hq' 

that "certain constants a^ and a 2 exist, and that £ k = 
v k (x) - aiv k (y) - a 2 (k = 1, K) conforms to a 

normal distribution of average 0 and variance 
(a 2 h(ai))". Here, a coefficient h(a^) is: 

20 h(ai) = 1 + ax 2 . . . (25) 

(Equation (25) is immediately derived from Equation 
(24) which gives the definitions of a^ and £ k , and a 
general property which a variance concerning a random 
variable has.) Besides, the value of the above 
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20 



variance o z can be estimated conveniently and 
satisfactorily precisely in practical use, by the 
method stated in, e.g., the foregoing section (III 
Case 1 of Applying Present invention to X-ray CT 
radiography; Dynamic CT radiography (To obtain 
Plurality of images of Reduced noises) ) . 

As understood from the above, if the constants a^ 
and a2 can be determined, it is possible to test the 
null hypothesis Hq' . Herein, it is actually sufficient 
that the optimum estimative values and a2~ of these 

constants are obtained. 

A known fitting method can be directly utilized 
for calculating such optimum estimative values of the 
constants a^ and a2 - In the ensuing description, 
therefore, a case of employing a linear least squares 
method shall be outlined as a typical concrete example 
of such a fitting method. In applying the linear least 
squares method to this embodiment, the sum-of-squares 
of in the above null hypothesis is merely defined as 
S (a) as follows : 



The value of the sum-of-squares S (a) depends upon a 
constant vector a = (a^, a2)/ that is, the values of 
the constants aj_ and a2 - When the constant vector a 
according to which the sum-of-squares S(a) takes the 
least value is calculated, the optimum estimative 
values ai~ and a2" in the sense of unbiased estimation 




(26) 



k = l 
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concerning the constants and a2 are obtained. 

Incidentally, various known methods can be utilized for 

the concrete computational method of the linear least 

squares method, and moreover, any of the known 

5 computational methods is very simple and requires a 

very slight computing time period. 

In this way, the optimum estimative values ai~, 

a 2 ~ of the constants a]_, a2 have been calculated, with 

the result that a residual defined by the following 

10 equation can be concretely computed: 

rjc~(X/ y) = v k (x) - a 1 "'v k (y) - a 2 " (27) 

Using the residual rj^, accordingly, the above null 

hypothesis Hq' can be said in other words as a null 

hypothesis Hq" being, in effect, equivalent thereto, 

15 that (x, y) (k = 1, • • * , K) conforms to a normal 

distribution of average 0 and variance (1 + (ai~) 2 )a 2 " . 

This is a concrete proposition which permits the 

% computations of testing to be actually executed. 

Further, the null hypothesis Hq' is said in other 

20 words by introducing vectorial expressions: 

a- = (ar,a 2 -> 1 _ ^ 

s = < r r' — * ztoj 

(where the vectors a and E, depend upon the combination 
x, y of the pixels) and employing a vector value 
function f defined by the following equation: 
25 f(a~, v(y)) = ai^v(y) + a 2 ~ ••• (29) 

Then, the null hypothesis Hq" becomes "v(x) = 
f(a~, v(y)) + £ (where £ conforms to a normal 
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distribution of average 0 and variance 

(1 + (ai^) 2 ) a 2 ) ", and this is in quite the same form as 
that of the null hypothesis Hq stated in the section 
(IX Most general form of Coherent filter in 
5 Embodiment) . It is accordingly apparent that this 

embodiment is one modification of the coherent filter 
according to the present invention. Incidentally, the 
function f(a~, v(y)) here signifies that parameter a 
representative of the partial volume effect which has 
10 been optimally adjusted for the pixel value v(y) of the 

pixel y and transformed so as to have the highest 
goodness of fit with the pixel value v(x) of the 
pixel x. 

Next, there will be explained a method in which, 
15 in this embodiment, a time-density curve is acquired by 

the coherent filter with the above null hypothesis Hq". 
Regarding a set R of pixels substantially corresponding 
\ to a part to-be-measured, computations to be stated 

below are executed for a certain pixel xeR included in 
20 the set R, relative to all pixels yeR being the 

elements of the set R. A residual (x, y) (k = 

1, •••/ K) is actually calculated using the above 
method. Subsequently, a critical probability p (x, y) 
or a weight w(p(x, y) ) in the case of abandoning the 
25 above null hypothesis Hq" that (x, y) (k = 1, 

K) conforms to a normal distribution of average 0 and 
variance (1 + (ai") 2 )a 2// is concretely computed in 
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accordance with Equation (6). Besides, a weighted 
average v'j^x) is computed in accordance with the 
following equation (15' ) so as to construct a time- 
density curve {<tfc, v'^fx) - v'i(x)> (k =1, 2, 
5 K) } at the pixel x. 

^ w (p (x, y)f ( a~,- v(y) ) 

V (x) = ^ = • • • (15' ) 

2^ w (p (x, y) ) 

yeR 

The time-density curve thus obtained is a 
measurement value which approximates {<tj c , A(d] c - d]_)>} 

10 (where A denotes an unknown coefficient) being the 

linear transformation of a true time-density curve 
{<tj cr d] c >} at the pixel x, and it has had random noises 
suppressed by the effect of the weighted average based 
on Equation (15). Besides, the pixel value vector of 

15 the other pixel y for use in the computation based on 

Equation (15) has had the influence of the partial 
* volume effect compensated as understood from the 

equation. Further, this embodiment has the property 
that "a temporal average is not used at all, and that a 

20 spatial average is computed using a weight based on the 

goodness of fit with the pixel x" as is the common 
feature of the coherent filter. Accordingly, the time- 
density curve which does not have a temporal resolution 
spoilt, which has the influence of the partial volume 

25 effect suppressed and which has the random noises 

suppressed can be obtained by this embodiment. 
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Incidentally, the scheme for acquiring the time-density 
curve in this way shall be especially termed "coherent 
regression method". 

Next, there will be concretely explained an 
5 example of the clinical utilization of a time-density 

curve in a dynamic image obtained in, e.g., dynamic CT 
radiography in X-ray CT for medical treatment. In this 
example of application, while a contrast medium is 
being rapidly injected into a blood vessel, the 

10 radiographing of dynamic CT or the like is performed, 

and the density change of the image of an artery 
existing in a human body tissue is measured as the 
time-density curve, thereby to intend the diagnosis of 
a blood flow movement in the tissue. 

15 In this example of application, the artery in the 

human body tissue is generally very fine in many cases, 
and hence, the image of the artery appearing on a 
tomogram based on the CT undergoes the partial volume 
effect. Further, random noises are, of course, 

20 contained in the image. With the prior-art method, 

therefore, it is difficult to obtain a satisfactorily 
precise time-density curve concerning the artery. When 
a time-density curve is daringly measured, it has been 
possible to obtain only measurement values <tk, 

25 (v^fx) - vi(x))> which somewhat approximate the linear 

transformation <tj c , AfD^ - D^)> of a true time-density 
curve <tj c , Dj c > concerning the artery (where denotes 



the pixel values (being scalar values) of a group of 
pixels corresponding to the image of the artery, at a 
time tfc, and k = 1, 2, K holds) . The measurement 

values contain random noises. Besides, the coefficient 
A remains unknown because of the influence of the 
partial volume effect . 

Therefore, when the above scheme according to the 
present invention is applied, it is possible to obtain 
measurement values <t] c , (v'] < (x) - v' i (x) ) > (k = 1, 
2, • • • , K) which satisfactorily approximate the linear 
transformation <t^, AfD^ - Di)>. On the other hand, 
veins which can be observed on the same tomogram 
include considerably thick ones. Regarding such a 
vein, accordingly, the satisfactorily good approximate 
values <tfc, (Jfc - Ji)> (k = 1, 2, K) of the time- 

density curve can be obtained by the prior-art method. 
Here, denotes the pixel values of a group of pixels 
corresponding to the image of the vein, at a time t^. 

Meanwhile, as to time-density curves concerning 
blood circulation, there has been known the inherent 
property of the proposition S: that "if a contrast 
medium density in blood is 0 (zero) at a time t]_, the 
areas under curves (AUCs) agree in the time-density 
curves <tj c , (d^ - di)> concerning any blood vessels d" . 
The "area under curve" mentioned here signifies the 
integral of the time-density curve <t) c/ (d^ - d^)> with 
respect to the time t. 
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Accordingly, the area-under-curve AUC (d) of the 
time-density curve <tj c , (d^ - d^)> concerning a certain 
blood vessel d can be approximately computed by, for 
example, the following equation: 

5 

AUC(d) - K f 1( (dk + 1 = dl) + (dk = dl) } (tk + 1 = tk) 

2 

k=l 

••• (30) 

Thus, an area-under-curve AUC (J) concerning the 
time-density curve {<t] c , (J^ - Ji)>} obtained as to the 
10 vein by the prior-art method can be computed using 

Equation (30) (J may be substituted into d) . Besides, 
regarding the artery, if a time-density curve {<tfc, 
(Dfc - Di)>} is known, an area-under-curve AUC(D) can be 
similarly computed using Equation (2 6) . Moreover, the 
15 following ought to hold in accordance with the above 

proposition S: 

AUC (D) = AUC (J) • • • (31) 

S In actuality, however, the time-density curve <t] c , 

(Dj£ - D]_)> is unknown, and hence, the area-under-curve 
20 AUC(D) cannot be computed. 

In contrast, the time-density curve <tj c/ (v'^tx) - 
v' i (x) ) > obtained by the scheme according to the 
present invention approximates <t] c , A(Dj c - D^)>, and 
the latter includes the unknown coefficient A. 
25 Therefore, an area-under-curve AUC (v' ) which can be 

concretely computed from the curve {<tj c , (v'^tx) - 
v']_(x))>} by using Equation (30) must be A times as 



t 
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large as the area AUC(D). That is: 

AUC(v') = A AUC(D) ••• (32) 

From Equations (31) and (32), the following 
relationship holds : 
5 A = AUC(v' ) /AUG (J) ••• (33) 

Since the right side of Equation (33) can be concretely 
computed using Equation (30) , the value of the 
coefficient A having been unknown can be concretely 
determined. Therefore, when a time-density curve <t) c , 
10 (v'] c (x) - v']_(x))/A> is constructed using the value of 

the coefficient A, it is nothing but approximating the 
time-density curve <t] c/ (D^ - Di)> of the artery. The 
method of constructing the time-density curve for which 
the value of the proportionality coefficient A having 
15 been unknown is determined, by using the area-under- 

curve in this manner, shall be called "AUC method" . 

As thus far described, in the clinical utilization 
% of a time-density curve in a dynamic image obtained by 

dynamic CT radiography or the like, the AUC method is 
20 further combined with the coherent regression method, 

whereby measurement values from which the influences 
of the partial volume effect and random noises are 
excluded and which do not include the unknown 
proportionality coefficient A are obtained concerning 
25 even the time-density curve of a fine artery the 

measurement of which has been difficult or impossible 
with the prior-art method. 
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Incidentally, the AUC method can, of course, be 
singly applied to the time-density curve <t^, (v'j c (x) - 
V i (x) ) > concerning the artery measured by the prior- 
art method, and the time-density curve for which the 
value of the proportionality coefficient A having been 
unknown is determined can be constructed ( though the 
influences of the random noises and partial volume 
effect cannot be excluded) . 

A time-density curve obtained by further combining 
the AUC method with the coherent regression method 
becomes, for example, as indicated by a sign (A) in 
FIG. 11 as to a certain pixel x. By the way, a time- 
density curve constructed by the prior-art method, to 
which the AUC method is singly applied, is also 
indicated at (B) in the figure. It is understood that, 
in the curve (B) , the influence of random noises is 
clearly seen, whereas in the curve (A) , the noises are 
sufficiently suppressed, with a temporal resolution not 
spoilt at all. 

Incidentally, although the scalar value d^ has 
been used as the density value of the time-density 
curve in the above description, this embodiment is not 
restricted to such a case. Needless to say, this 
embodiment can be easily expanded and applied also to a 
case, for example, where individual images constituting, 
a dynamic image are color images, more generally the 
combination of a multiple sorts of images, and where 
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the density value of a time-density curve is expressed 
as a vector value . 

(X Supplementary items of Embodiments) 

Now, the supplementary items of the embodiments 
5 thus far described will be described* 

(X-l General supplementary items) 

First, in the foregoing embodiments, the "critical 
probability" p (x, y) in the case of abandoning the null 
hypothesis H has been supposed as means for quantifying 

10 the "goodness of fit" termed in the present invention, 

but the present invention is not restricted to such an 
aspect. As stated at the beginning, the "goodness of 
fit" may be a numerical "index" which represents 
whether or not the pixels x and y are similar in any 

15 sense. 

Besides, in this regard, the x~square test method 
has been employed above in order to calculate the 
critical probabilities in the case of abandoning the 
several null hypotheses (refer to Equations (6) and 

20 (9)), but the present invention is not restricted 

thereto. In other words, the present invention is not 
restricted to only the aspect of utilizing the 
"statistical testing method" or the "x~square test 
method" which is one sort of the concrete aspects 

25 thereof. 

Further, in the foregoing embodiments, the weight 
functions totaling two sorts, wl (Equation (10) into 
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which Equation (6) or (9) has been substituted) and w3 
(Equation (14)) have been mentioned as the concrete 
form of the weight function w(p(x, y) ) , but the present 
invention is not restricted thereto. As already 
stated, the property which is to be satisfied by the 
weight function w being the function of the critical 
probability p(x, y)e[0, 1] is only the fact that the 
weight function w(t) is the non-negative monotonously 
increasing function defined in a domain te[0, 1]. 
Besides, most generally speaking, the weight function w 
may be the "non-negative monotonously increasing 
function which concerns the goodness of fit". 
(X-2 More general aspect concerning Goodness of fit) 

The "goodness of fit" termed in the present 
invention is not restricted to the scheme as stated 
above wherein it is numerically given in accordance 
with the critical probability p(x, y) in the case of 
abandoning the null hypothesis. Accordingly, the 
"weight" is not restricted, either, to the scheme 
wherein it is calculated by exerting the weight 
function w on the critical probability p(x, y) . In the 
ensuing description, more general aspects concerning 
such goodness of fit, etc. shall be explained along 
practicable examples. 

First, to sum up generally, the present invention 
adopts a construction wherein, regarding a certain 
proposition appropriately set in conformity with the 



purpose of image processing, a certainty at which the 
proposition holds in an image to-be-processed is 
numerically given by an appropriate scheme conforming 
to the purpose, and the numerical value is subsequently 
applied to the image processing. By way of example, in 
order to numerically give the goodness of fit between a 
pair of pixels x and y, a more moderate scale may well 
be constructed in the present invention. As a 
practical example thereof, there is considered a 
construction wherein the goodness of fit is numerically 
given by utilizing a membership function in fuzzy logic 
or the like. 

(X-2-1 Example of Applying Present invention to 
Pattern recognition) 

Here, as a practicable example a more general 
aspect concerning the goodness of fit, etc., there will 
be explained one example in the case of applying the 
present invention to so-called "pattern recognition" 
which identifies and extracts a featuring pattern on an 
image . 

In this case, on the basis of the feature of a 
pattern which is to be identified, a function m(x) 
which represents the certainty of the proposition that 
"a pixel x is a pixel constituting the pattern" is 
defined beforehand. Here, m(x)e[0, 1] is held. In 
more detail, the function m(x) is a function which 
calculates a numerical value representative of the 
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certainty of the proposition from a vector value v(x) 
being the pixel value of the pixel x, the pixel value 
v(y) of each pixel y included in a set N(x) of pixels 
lying in the vicinity of the pixel x, a vector x itself 
representative of the position of the vector, and so 
forth. Accordingly, the function m(x) needs to be 
appropriately set in accordance with the difference 
between the pattern which is to be extracted and a 
pattern which is not to be extracted. Besides, a 
weight function is defined by: 



Besides, when a set X of pixels as to which the 
function m(x) is larger than a certain "threshold 
value" is constructed, the region of the pixel x 
corresponding to the pattern is obtained. 

Now, when an image P to be processed is given, 
w(T, m(x)) in which a concrete threshold value T and 
the function m(x) are substituted into the above weight 
function is computed as to all the pixels x of the 
image. (Thus, the set X whose membership function is 
w(T, m(x)) (that is, the set of the pixels x as to 
which w(T, m(x) ) = 1 holds) is defined, and this set X 
is nothing but the region occupied by the pattern on 
the image P.) Further, a new image P' is generated as 
stated below. The pixel x of the image P' has v(x) as 
its pixel value when the function w(T, m(x)) = 1 holds, 




(34) 
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and it has 0 (zero vector) as its pixel value in any 
other case. The image P' thus created is such that the 
pixel values of pixels which do not correspond to the 
pattern to-be-identified are erased from the image P 
5 by 0. 

Of course, the processing steps of computing the 
function m(x), computing the weight w(T, m(x)) and 
constructing the image P' may well be respectively 
executed in sequence, or they can also be executed by 
10 software in which these processing steps are 

integrated. 

In this example, the function m(x) numerically 
gives the ''goodness of fit" between the pattern and the 
pixel x by checking the feature of the pattern, the 
15 vector value v(x) being the pixel value of the pixel x, 

and the value of the vector x itself. Besides, the 
weight w(T, m(x)) is calculated by the weight function 
S defined by Equation (34), and the new image P' being 

the processed result is constructed using this weight 
2 0 and the image P. 

As a further practicable example concerning this 
processing, there will be described the construction 
of image processing for extracting a fine line in an 
aerial photograph, the microphotograph of a metal 
25 surface, or the like. 

First, there will be explained a method for 
constructing an index which numerically represents the 



certainty of the proposition that "a pixel x is a pixel 
constituting the fine line". 

A set N(x) of pixels lying in the vicinity of the 
pixel x is defined beforehand. Typically, the set N(x) 
may be defined as a rectangular area which centers 
around the pixel x. Also defined is a multidimensional 
vector value v(x) in which the pixel values (scalar 
values or vector values) of the individual pixels in 
the set N(x) are arrayed. The number of dimensions of 
the vector is denoted by K. Typically, the number K is 
preferably set at about 25 - 169. 

Subsequently, a normalization function nrm being a 
K-dimensional vector function which is auxiliarily used 
in later computations is defined as stated below. The 
normalization function nrm is given for any K- 
dimensional vector v subject to v ^ 0, as follows: 



where k = 1, 2, K holds. Owing to this 

normalization function, it is guaranteed that the 
following holds for any K-dimensional vector v (as long 
as the vector v is not a zero vector) : 




(35) 



K 

2] nrmk(v) = 0 
i = l 



Z e rm : 

i = l 




Subsequently, one pixel x is determined at will, 
and typical examples of images each of which has a 
background of 0 (zero) and includes only one fine 
straight line passing through the pixel x are collected 
in a suitable number J and have image Nos. 1, 2, J 
assigned thereto. The value of v(x) in the image of 
image No. j is computed, thereby to concretely 
construct a vector r(j) in accordance with the 
following equation: 

ri ( j ) = nmr j (v (x) ) 
Such computations are executed as to j = 1, 2, •*•, J. 
(The J vectors r(j) thus constructed shall be called 
"pattern vectors' 7 below. Incidentally, the pattern 
vector corresponds to the "pixel value separately 
constructed" termed in the present invention.) 

Besides, a function m(x) is defined by the 
following equation: 

m(x) = max{r ( j ) •nrm(v(x) ) I j = 1, 2, • • •, J} 
Here, "•" denotes the inner product between the K- 
dimensional vectors, and "max" denotes a function for 
deriving the maximum value in the elements of the set. 
The function m(x) thus constructed takes a large value 
in a case where the pixel x matches with or is similar 
to any of the pattern vectors r(j) (j = 1, 2, • • • , J) , 
and it takes a small value in any other case. 
Therefore, this function m(x) is the index which 
numerically represents the certainty of the proposition 
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that "the pixel x is a pixel constituting the fine 
line", that is, a function which represents the 
"goodness of fit". Further, a weight function is 
defined by: 



Here, T denotes a constant which serves as a threshold 
value, and as which a suitable value is set. 

Now, when an image P to be processed is given, 
w(T, m(x)) in which a concrete threshold value T and 
the function m(x) are substituted into the above weight 
function is computed as to all the pixels x of the 
image. (Thus, the set X whose membership function is 
w(T, m(x)) (that is, the set of the pixels x as to 
which w(T, m(x) ) = 1 holds) is defined, and this set X 
is nothing but a region occupied by the pattern on the 
image P . ) 

Further, a new image P' is generated as stated 
below. The pixel x of the image P' has v(x) as its 
pixel value when the function w(T, m(x)) = 1 holds, and 
it has 0 (zero vector) as its pixel value in any other 
case. The image P' thus created is such that the pixel 
values of pixels which do not correspond to the pattern 
to-be-identified are erased from the image P by 0. 

Besides, w(T, t) may well be constructed as a 
fuzzy discrimination function, for example: 




(34 mentioned again) 



w(T, t) = 4 
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(3(1)2 _ 2(-) 3 )C for 0 < - < 1 
T T T 

0 for — < 0 

T 

1 for - > 1 

T 



(where C denotes a positive constant), and the pixel x 
of the image P' may well be constructed so as to have 
w(T, m(x)V(x)) as its pixel value. This lowers the 
contrast of a pixel which is judged as not considerably 
corresponding to the "fine line", by the fuzzy logic, 
thereby to produce the effect of constructing the image 
P' so that the part of the "fine line" may stand out in 
relief. 

(X-3 Supplementary items concerning Estimation method 
for Variance a 2 ) 

Next, supplementary description will be made on a 
method of estimating the variance a 2 (or standard 
deviation a) for use in the above equation (6) or (14) . 

In the foregoing "Case of Applying Present 
invention to X-ray CT radiography", in the light of 
the background that the expected value E[a 2 ] of the 
variance obtained by the above equations (6) and (7) is 
proper for all the pixels x on all the K images, this 
expected value E[a 2 ] has been indiscriminately used for 
all the pixels x. However, there is also a case where 
such a premise is not satisfied, that is, a case where, 
as to the respective pixels x (which are x^, and 
x F here), peculiar variances o 2 (xi) , and a 2 (x F ) 

are generally supposed (even if the shapes of 
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distributions are identical (for example, a Gaussian 
distribution) ) • 

On such an occasion, the expected values 
E[a 2 (xi)]/ ••*/ and E[a 2 (xp)] of variances may well be 
5 calculated individually for the respective pixels x in 

accordance with the above equations (6) and (7). 
Thereafter, the expected values E[a 2 (xi)], and 
E[a 2 (xp)] individually calculated in this manner may be 
substituted into Equation (3) in correspondence with 
10 the respective cases (= in correspondence with the 

pixels xi, •••/ and xp noticed) so as to obtain v' (x) = 
(v'!(x), v' 2 (x), v' K (x)). 

Besides, in a case where noises contained in the 
pixel values of all pixels are supposed to have a 
15 common variance a 2 , a scheme as stated below can be 

adopted for estimating the variance a 2 . First, as to 
the set Q of all pixels x (= {x^, xp}), the 

) average (a 2 ) 4 " of the expected values E[a 2 (x)] of 

variances as calculated by the above equations (7) and 
20 (8) is obtained, and it is used as the first 

approximation of the estimative value of a 2 in Equation 
(6) . That is: 

( q2 > + = |H 2>[a 2 (x)] ' (35) 

|n| xen 

25 In the above equation, |Q| denotes the number of the 

elements of the set Q (here, |Q| = F holds) . 
Subsequently, M is let denote a set of pixels x as to 
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which the expected values E[a^(x)] lie within, for 
example, several times the above average value (a 2 ) + . 



the pixels x which are the elements of the set M, it 
becomes : 



This average can be utilized as the more likely 
estimative value of a. 

Such processing is effective in, for example, a 
case as concretely stated below, that is, a case where 
two static images Gl and G2 constituting a dynamic 
image exist as shown in FIG. 12(a), and where the 
difference between the images Gl and G2 is taken as 
shown in FIG. 12 (b) , images Zl and Z2 somewhat deviate 
from each other (in other words, the image Zl has moved 
as the image Z2) . 

The variance (a 2 ) + of the above equation (35) 
corresponds to a variance calculated on the basis of 
all the pixels x (= Q) constituting an image G3 as 
shown in FIG. 12 (b) . In this case, however, two peaks 
Cl' and C2' / for example, ascribable to parts £1 and C2 
where the images Gl and G2 are not superposed on 
each other are caused to appear in a probability 
distribution diagram concerning the image G3 as shown 
in FIG. 12(c). Therefore, the variance (a 2 ) + is 
overestimated. It is accordingly inappropriate to use 



When (a 2 ) ++ 



is let denote a new average taken for only 




(36) 



the variance (a 2 ) + left intact, for all the static 
images. On the other hand, as to a part Cm shown in 
FIG . 12(b), a central peak C^ r in FIG . 12(c) is 
considered to correspond thereto. It is therefore more 
preferable to calculate a variance concerning a 
probability distribution diagram being in a state where 
the peak £m' is left by removing the two peaks £l' and 
£2' • Such a variance really corresponds to the (a 2 ) ++ 
in the above equation (32) . 

In other words, the above processing that the 
average is taken for only the set M of "pixels x as to 
which the expected values E[a 2 (x)] lie within several 
times the average (a 2 )" 1 "" is nothing but processing in 
which an average is taken by extracting only the peak 
0m' in FIG. 12(c). As a result, the calculation of the 
average (a 2 ) ++ is equivalent to that of the variance 
concerning the probability distribution diagram in the 
state where the peaks £l' and £2' have been removed in 
FIG. 12 (c) . 

In this case, accordingly, the average (a 2 ) ++ is 
substituted into the above equation (5), etc., whereby 
the more likely values of the transformed pixel values 
v' (x) = (v'i(x), v'2< x )/ v' K (x)) concerning the 

respective pixels x^, and xp can be obtained. 

Incidentally, regarding the estimation of noises, 
in general, a method is also effective wherein a 
subject of uniform brightness is put at a position 
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which is to lie, e.g., at the corner of an image, and 
the distribution of the values of a set of pixels at 
the corresponding part is measured, thereby to estimate 
the distribution of the noises. 
5 In the present invention, the method of estimating 

the distribution or variance of noises may fundamen- 
tally be based on any scheme, and it is not restricted 
to several schemes stated here. In practice, a scheme 
is to be properly adopted wherein the distribution or 

10 variance can be most preferably calculated using 

empirical information, an imaging process theory, 
measurement values or actually measured values, 
etc. available in each of various embodiments. 
Incidentally, on this occasion, it is favorable for 

15 simplifying the construction of an apparatus and 

enhancing a processing speed that the estimation method 
is constructed as easily as possible, in conformity 
with a practical precision required in each embodiment. 
(X-4 Further embodiments of Present invention) 

20 Lastly, supplementary description will be made on 

further embodiments of the present invention as have 
not been referred to in the foregoing. 
(X-4-1 Embodiment as Preprocessing of Data 
compression) 

2 5 In general, when an image Q' in which random 

noises are added to a certain image Q is created, the 
image Q' has a much larger information content as 
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compared with the image Q (an exception is a case where 
the original image Q contains very large noises) . 
Accordingly, if such an image Q' is to be communicated 
or to be recorded on a recording medium, an information 
5 content for describing the noises having no substantial 

significance uses the number of bits in the communi- 
cations or recording wastefully. Conversely speaking, 
when an image Q' ' with the noises suppressed is created 
from the image Q' and is handled for the communications 

10 or recording, an image having the same significance 

(and being vivid) can be expressed with a small 
information content . 

In general, in case of handling an image in 
communications or recording, the image is often 

15 subjected to so-called "compression processing". This 

processing is intended to reduce an information content 
to-be-transmitted or -recorded to the utmost, thereby 
* to realize high-speed communications or large-quantity 

recording (for one recording medium) . In practical 

20 use, it is in the case of a "dynamic image" that the 

information content is to be especially reduced. These 
are facts which are extensively known in irreversible 
image compression technology. 

In view of these items, the number of bits (the 

25 information content) necessary for the communications 

or recording can be made smaller by compressing the 
image Q' ' with the noises suppressed, rather than by 
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compressing the image Q' containing the noises. 

Therefore, when the image Q' to be communicated or 
recorded is subjected to preprocessing for suppressing 
the noises by applying the coherent filter according to 
5 the present invention (concretely, a "preprocessing 

apparatus" is constructed by hardware or software, and 
the same applies below) , the image Q' ' with the noises 
suppressed is obtained without spoiling a temporal 
resolution or a spatial resolution, and this image Q' ' 

10 can be compressed at a high compression rate, whereby a 

time period and a memory capacity necessary for the 
communications or recording can be sharply saved, 
(X-4-2 Embodiment as Preprocessing of Two-dimensional 
or three-dimensional regional extraction processing) 

15 Processing for extracting specified regions from 

the two-dimensional image of an aerial photograph or 
the like or from the three-dimensional distribution 
\ image based on the MR I equipment or the like, is often 

executed. By way of example, the regions of only roads 

20 are extracted regarding the former aerial photograph, 

while the regions of only blood vessels are extracted 
regarding the latter three-dimensional distribution 
image of a human head as based on the MRI equipment or 
the like, and the extracted regions are subjected to 

25 rendering so as to display them as computer graphics. 

Especially, the blood vessels in the latter can be 
relatively rotated or observed in various directions by 
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way of example, and the radiography thereof is 
recognized as important technology in the diagnoses of 
a cerebrovascular disease, etc. 

In case of such regional extraction, the simplest 
5 and most fundamental processing is so-called "threshold 

processing". The threshold processing is processing in 
which individual pixels x constituting the two- 
dimensional image or three-dimensional distribution 
image, or the like are classified, depending upon if 

10 their scalar values are larger than a certain 

predetermined threshold value. Besides, when a set of 
pixels consisting of only the pixels x which fall 
within a specified class is subjected to rendering, an 
image having undergone the regional extraction, for 

15 example, a "three-dimensional image consisting of only 

blood vessels'" can be displayed (incidentally, various 
schemes other than the "threshold processing" have also 
been used as prior-art techniques) . 

Meanwhile, in such regional extraction processing, 

20 it is desirable that noises are little contained in the 

original image (two-dimensional image or three- 
dimensional distribution image mentioned above, or the 
like) . This is because the state occurs that, when 
noises exist in the image, a pixel which lies in a 

25 place which ought not to be, for example, a blood 

vessel has a pixel value as if it were the blood 
vessel, whereas a pixel which lies in a place which 



ought to be blood vessel has a pixel value as if it 
were not the blood vessel- In this state, accordingly, 
the step of removing such a pixel erroneously 
classified needs to be performed by complicated 
postprocessing. Besides, in the actual situation, the 
"smoothing processing" as stated in the section of the 
prior art is generally performed in order to eliminate 
the noises on the original image. In this case, 
however, there is the problem that a spatial resolution 
is spoilt, so the fine structure of the region to-be- 
extracted (including, for example, fine blood vessels) 
are lost. 

Therefore, when the original image, namely, the 
two-dimensional image or three-dimensional distribution 
image, is subjected to preprocessing for suppressing 
the noises by applying the coherent filter according to 
the present invention, the noises of the image can be 
suppressed without spoiling the spatial resolution, 
with the result that precise regional extraction 
processing can be carried out. 

Incidentally, the embodiment similar to such 
extraction processing has already been described in the 
foregoing section "Case of Constructing Pixel value 
v(x) from Single image" (refer to FIG. 8(c)). 
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CLAIMS 

1. An image processing apparatus characterized by 
comprising: 

decision means for deciding a similarity between a 
5 first pixel and a second pixel which constitute an 

image, by testing; and 

average means for subjecting the first pixel and 
the second pixel to weighted averaging on the basis of 
a decided result by the decision means. 
10 2 . An image processing apparatus characterized by 

comprising: 

numerical means for numerically giving a weight 
which is determined by a similarity between a first 
pixel and a second pixel constituting an image; and 
15 average means for subjecting values of the first 

pixel and the second pixel to weighted averaging by 
using the weight numerically given. 

3. An image processing apparatus as defined in 
claim 1 or 2, characterized in that the average means 

20 obtains a new pixel value concerning the first pixel, 

as a weighted average. 

4. An image processing apparatus as defined in 
claim 1, characterized in that the average means 
includes determination means for determining the 

25 weighting factor on the basis of the decided result, 

and multiplication means for multiplying values of the 
first pixel and the second pixel by the weighting 



factor . 

5. An image processing apparatus as defined in 
any of claims 1, 2, 3 and 4, characterized in: 

that the average means multiplies the second pixel 
value by a large weighting factor by the multiplication 
means when the similarity between the first pixel and 
the second pixel is high, and multiplies the second 
pixel value by a small weighting factor when the 
similarity is low; and 

that it subjects the weighted first pixel value 
and second pixel value to the weighted averaging. 

6. An image processing apparatus as defined in 
claim 1, characterized in that the decision means 
decides the similarity between the first pixel and the 
second pixel on the basis of a plurality of images 
obtained in such a way that an identical subject is 
photographed at different times by employing an imaging 
equipment . 

7. An image processing apparatus as defined in 
claim 1 or 6, characterized in that the decision means 
decides the similarity between the first pixel and the 
second pixel on the basis of a plurality of images 
obtained in such a way that an identical subject is 
photographed at different photographing conditions by 
employing an imaging equipment. 

8. An image processing apparatus as defined in 
any of claims 1, 6 and 7, characterized in that the 



decision means decides the similarity between the first 
pixel and the second pixel on the basis of a plurality 
of images obtained by photographing an identical 
subject with an imaging equipment, and under different 
processing conditions . 

9. An image processing apparatus as defined in 
claim 1, characterized in that the decision means 
decides the similarity between the first pixel and the 
second pixel by using a result of a comparison which is 
made between a vector value constructed by arraying 
scalar values of the first pixel in a single image 
obtained by photographing a subject with an imaging 
equipment, and a vector value constructed by arraying 
scalar values of the second pixel. 

10. An image processing apparatus as defined in 
any of claims 1, 6, 7 and 8, characterized in that, 
using a result of a comparison which is made between a 
first vector constructed by arraying corresponding 
pixel values of the plurality of images in coordinates 
of the first pixel, and a second vector constructed by 
arraying corresponding pixel values of the plurality of 
images in coordinates of the second pixel, the decision 
means decides the similarity between the first pixel 
and the second pixel in a predetermined one of the 
plurality of images. 

11. An image processing apparatus as defined in 
claim 1 or 2, characterized by further comprising image 



average means for averaging a first image obtained by 
subjecting the image to averaging processing with the 
average means, and a second image obtained by 
subjecting the image to processing different from the 
averaging processing. 

12. An image processing apparatus as defined in 
claim 1 or 2, characterized by further comprising image 
average means for averaging a first image subjected to 
averaging processing by the average means, and a second 
image before subjecting the first image to the 
averaging processing . 

13. An image processing apparatus characterized by 
comprising: 

determination means for determining and 
numerically giving a similarity between a first pixel 
and a second pixel which constitute an image, by 
statistical testing; and 

average means for subjecting values of the first 
pixel and the second pixel to weighted averaging on the 
basis of the numerical value. 

14. An image processing apparatus characterized by 
comprising : 

numerical means for numerically giving a 
similarity between a first pixel and a second pixel 
which constitute an image, by statistical testing; and 

average means for averaging values of the first 
pixel and the second pixel when the numerical 
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similarity is high in the numerical means, and for not 
averaging the first pixel value and the second pixel 
value when the determined similarity is low. 

15. An image processing method characterized in: 

5 that vector values each of which is constructed by 

arraying scalar values are calculated as to first and 
second pixels which constitute an image; and 

that, in quantifying a similarity between a value 
of the first pixel and a value of the second pixel 

10 constructed separately from the first pixel value and 

then constructing a new pixel value for the first pixel 
by utilizing the second pixel value, the new pixel 
value is constructed by making contribution of the 
second pixel large when the similarity is high, and by 

15 making contribution of the second pixel small when the 

similarity is low. 

16. An image processing method characterized in: 

^ that a first pixel value and a second pixel value 

each of which is a vector value constructed by arraying 
20 scalar values are constructed as to a first pixel and a 

second pixel which constitute an image; 

that a similarity between the first pixel value 
and the second pixel value is quantified; and 

that, in constructing a new pixel value for the 
25 first pixel by utilizing the second pixel value, 

the new pixel value is constructed by making 
contribution of the second pixel large when the 
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similarity is high, and by making contribution of the 
second pixel small when the similarity is low. 

17. An image processing method characterized in: 
that a first pixel value and a second pixel value 

5 each of which is a vector value constructed by arraying 

scalar values are constructed as to a first pixel and a 
second pixel which constitute an image; 

that a similarity between the first pixel value 
and the second pixel value is quantified; and 

10 that, in constructing a new pixel value for the 

first pixel, by exerting a weight function which is a 
function of the similarity, on the similarities 
calculated for the individual second pixel values, 
thereby to determine weights of the individual other 

15 pixel values, and by subsequently calculating a 

weighted average of the other pixel values on the basis 
of the weights, 

the new pixel value is constructed by making the 
weight large so as to make contribution of the other 

20 pixel large in the weighted average, when the 

similarity is high, and by making the weight small so 
as to make contribution of the other pixel small in the 
weighted average, when the similarity is low. 

18. An image processing method as defined in 

25 claim 17, characterized in that the weight function is 

a non-negative monotonously increasing function which 
concerns the similarity. 
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19. An image processing method as defined in 
claim 17, characterized in: 

that, when the first pixel and the second pixel 
are respectively denoted by x and y, the first pixel 
5 value and the second pixel value are respectively 

denoted by v(x) = (v x (x) , v 2 (x) , v K (x)) and v(y) = 

(vi (y) , V2 (y) t v K (y)), and the new pixel value to 

be constructed is denoted by v' (x) = (v' i (x) , 
v'2(x), ••*, v' K (x)), and when p(x, y) denotes the 
10 similarity, w denotes the weight function, and 

w(p(x, y) ) denotes the weight obtained by exerting the 
weight function w on the similarity, processing for 
taking the weighted average is expressed by: 

X v k (y)w(p(x, y) ) 

15 v' k (x, = y6N W - - 

2^ w(p(x, y) ) 

yeN(x) 

(where N(x) denotes a range in which the other pixel is 
included) . 

20. An image processing method as defined in any 
of claims 15 through 17, characterized in that noises 

20 of the image are suppressed by constructing the new 

pixel value. 

21. An image processing method as defined in any 
of claims 15 through 17, characterized in that pattern 
recognition in the image is performed by constructing 

25 the new pixel value. 

22. An image processing method as defined in any 
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of claims 15 through 17, characterized in that the 
similarity is quantified on the basis of a critical 
probability which is obtained by applying a statistical 
testing method to scalar values respectively 
constituting the first pixel value and the second pixel 
value . 

23. An image processing method as defined in 
claim 16 or 17, characterized in that the second pixel 
is selected from within a predetermined region around 
the first pixel. 

24. An image processing method as defined in 
claim 23, characterized in that the predetermined 
region is a rectangular area which centers around the 
first pixel. 

25. An image processing method as defined in 
claim 16 or 17, characterized in that the image 
consists of a plurality of images, and that the pixel 
values are vector values which are obtained by arraying 
scalar values of pixels of the same points throughout 
the plurality of images. 

26. An image processing method as defined in 
claim 25, characterized in that the plurality of images 
are a plurality of static images which constitute a 
dynamic image. 

27. An image processing method as defined in 
claim 2 6, characterized in that the dynamic image is a 
dynamic CT image which is acquired by a medical image 



diagnostic equipment . 

28. An image processing method as defined in 
claim 16 or 17, characterized in that the image is a 
plurality of sorts of images concerning an identical 
subject, and that the pixel values are vector values 
which are obtained by arraying scalar values of pixels 
of the same points throughout the plurality of sorts of 
images . 

29. An image processing method as defined in 
claim 28, characterized in that the image is an image 
which is acquired by multiwindow radiography in a SPECT 
equipment or PET equipment, and that the plurality of 
sorts of images are respectively different images which 
are based on gamma rays respectively emitted from 
different radioisotopes . 

30. An image processing method as defined in 
claim 28, characterized in that the image is a color 
image, and that the plurality of sorts of images are 
respectively different images which are acquired by 
decomposing light into three primary colors. 

31. An image processing method as defined in any 
of claims 15 through 17, characterized in that the 
image is a single image, and that the pixel values are 
constructed from the single image. 

32. An image processing method as defined in 
claim 31, characterized in that the pixel values are 
vector values which are obtained by arraying a scalar 
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value of a certain pixel on the single image and scalar 
values of pixels lying in the vicinity of the certain 
pixel • 
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ABSTRACT 
This invention consists in an image processing 
method in which the similarity between a first pixel 
and a second pixel constituting an image is numerically 
5 given by statistical testing, and the values of the 

first and second pixels are averaged when the 
similarity numerically given is high, whereas the first 
and second pixel values are not averaged when the 
determined similarity is low, whereby an image with 
10 noises suppressed is obtained without spoiling a 

spatial resolution as well as a temporal resolution. 
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